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ABSTRACT 


A  theoretical  model  has  been  developed  to  predict  acoustic 
radiation  from  a  submerged,  finite-cylindrical  shell  internally 
stiffened  by  ring  frames  and  with  flat,  movable  end  caps.  Iterative 
techniques  are  developed  to  solve  the  associated  fluid-structure 
coupled  systems  over  a  moderate  range  of  frequencies  which  encompasses 
several  fluid-loaded  resonances.  We  demonstrate  that  such  techniques 
generate  essentially  the  exact  solution  in  a  reasonably  small  number 
of  iterations.  Comparison  of  computational  costs  to  those  for  the 
exact  result  shows  iteration  to  be  economically  advantageous  for  all 
frequencies  considered. 
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I.  INTRODUCTION:  THE  MODEL  IN  CONTEXT 

A.  Statement  of  the  Problem 

The  prediction  of  acoustic  radiation  from  a  general  submerged 
body  requires,  as  an  intermediate  step,  that  two  coupled  systems  of  equa¬ 
tions  be  solved  simultaneously  for  the  surface  response  and  radiation 
loading.  One  of  the  systems  describes  the  in-vacuo  dynamics  of  the 
structure  and  so  appears  driven  by  unknown  fluid  pressures  as  well  as 
by  known  applied  forces;  the  second  describes  the  body-surface  structure¬ 
less  dynamics  of  the  fluid  and  in  exact  formulations  is  usually  provided 
by  the  Helmholtz  integral,  which  also  later  yields  far-field  pressures 
once  its  integrand  has  been  fully  determined.  While  the  generation  of  the 
latter  fluid  system  may  sometimes  be  expected  to  dominate  computational  costs, 
a  non-negligible  part  of  the  total  expense  must  always  be  earmarked  to  the 
coupled  systems'  solution  process  itself.  The  fluid  system  is  again 
mostly  to  blame  since  its  effective  contribution  to  the  total  impedance 
is  a  fully  populated  matrix  in  contrast  to  that  of  the  in-vacuo  structure, 
which  is  typically  heavily  banded  and  therefore  relatively  easy  to  invert. 

A  major  objective  here  is  the  development  and  implementation  of  cost- 
effective  iterative  methods  to  solve  fluid- structure  systems  which  arise 
in  the  analysis  of  a  finite,  frame-stiffened  cylindrical  shell  with  movable 
end  caps.  The  project  is  in  part  a  follow-up  to  that  of  Ref.  1,  where  we 
established  that  similar  iterative  techniques  worked  for  the  force-driven 
spherical  shell.  The  spherical  structure,  however,  is  a  separable  geometry 
which  besides  is  modally  compatible  with  its  external  medium.  The  iterations 
performed  where  thus  applied  to  scalar  quantities  and,  though  generally 
successful,  left  unanswered  the  fundamental  question  regarding  feasibility 
of  application  to  general  nonseparable  cases  with  fully  populated  matrices. 

It  will  be  recalled  that  when  applicable,  iteration  methods  for  the  solution 
of  a  single  linear  system  traditionally  involve  an  order  of  magnitude  fewer 
operations  than  do  exact  methods. 
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A  second  objective  is  to  apply  the  model  to  an  actual  structure 
and  to  compare  exact  far-field  predictions  to  existing  measurements. 

Summarizing  the  work,  Chapter  II  contains  the  development  of  modal 
structural  and  fluid  systems  for  a  submerged  finite  cylinder,  and  Chapter  III 
that  of  the  several  iteration  algorithms  intended  to  solve  them.  In  Chapter  IV 
the  theory  is  applied  to  an  actual  situation  of  interest  and  exact  and  iterated 
solutions  are  discussed  and  compared  to  measurement. 

B.  Background  ' 

The  following  is  a  brief  review  of  recent  developments  in  the  area 
of  vibration  of  cylindrical  structures,  framed  or  unframed,  submerged  in  a 
heavy  fluid  which  at  times  may  substantially  change  the  character  of  the  in- 
vacuo  response.  The  use  of  iteration  in  past  studies  is  also  discussed  briefly. 

B. 1  Treatment  of  the  External  Fluid  Space 
B.l.a  "Exact"  Models 

As  earlier  stated,  rigorous  treatments  of  the 
external  fluid  have  tended  to  rely  on  the  Helmholtz  integral,  the  alternative 
being  a  finite-element  or  finite-difference  partitioning  of  a  chosen  fluid 
volume  enclosing  the  body  and  the  enforcement  of  a  radiation  condition  on  the 
outmost  elements.  The  Helmholtz  integral  has  the  advantage  that  it  yields 
an  integral  equation  which  requires  pressure  and  displacements  values  only 
on  the  body  surface  —  the  values  which  completely  determine  an  interaction 
problem. 


A  major  disadvantage  however,  is  that  by  itself  the  Helmholtz  integral 

cannot  be  solved  at  the  eigenfrequencies  of  its  kernel,  so  that  fluid  im- 

2 

pedances  are  not  always  immediately  available  and  therefore  neither  are  values 
of  the  fluid  pressure  which  would  accompany  a  specified  boundary  motion. 

This  disadvantage  is  only  apparent, however,  since  the  fluid  eigenproblem 
is  removed  upon  coupling  to  an  interior  problem  posed  for  the  body  cavity 
which  complements  the  entire  medium,  and  techniques  such  as  those  discussed 
in  Refs.  3,4  may  be  employed  for  a  general  impedance  calculation;  the 
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i 
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methods  of  Refs.  5  and  6  achieve  the  same  end  by  coupling  the  Helmholtz 
integral  to  the  related  equation  resulting  from  differentiation.  More¬ 
over,  in  an  interaction  problem  this  apparent  disadvantage  is  entirely 

7 

fictitious  since,  as  pointed  out  by  Huang  ,  fluid  impedances  need  then 
never  be  calculated  as  a  bridging  step.  In  effect,  the  fluid  is  then 
automatically  coupled  to  an  "interior"  structural  system  and  the  Helmholtz 
eigenproblem  does  not  come  up. 

Still,  a  number  of  useful  and  interesting  studies  have  been  ex¬ 
clusively  devoted  to  the  calculation  of  fluid  impedances  for  the  cylin¬ 
drical  geometry,  and  thus  to  the  mapping  out  of  the  radiation  loading 

8 

over  such  structures.  Chertock  for  example,  has  developed  a  general 

computational  scheme  for  obtaining  solutions  of  the  Helmholtz  integral 

equation  corresponding  to  given  or  assumed  velocity  distributions.  The 

method  is  iterative  and  he  observes  that, while  it  converges  easily  for 

low  frequencies,  higher  frequencies  demand  certain  modifications,  e.g. , 

averaging  of  two  previous  iterates  for  input  in  the  current  step.  Slow 

or  nonconvergence  of  iterative  techniques  applied  just  to  the  fluid  system 

has  been  connected  in  Ref.  1  to  the  Helmholtz  eigenfrequencies  by  recalling 

the  convergence  constraints  on  the  standard  Neumann  series  for  a  Fredholm 
9 

equation  . 

Sandman  has  investigated  the  finite  but  baffled  cylinder'*'0,  and  the 
finite  unbaffled  cylinder  with  zero  prescribed  motion  at  the  two  flat  ends. 
His  objective  was  to  determine  modal  radiation  loading  magnitudes  corre¬ 
sponding  to  prescribed  modal  motions,  with  results  which  could  be  layed 
out  over  those  for  the  classical  infinite  cylinder*  solution  to  determine 
the  effect  of  the  finite  ends.12  Similarly,  Kozyrev  and  Shenderov13 ' 14 
have  calculated  fluid  impedances  for  a  finite  cylinder  with  rigid  and  soft 

3 

caps  using  a  finite-difference  approach;  they  applied  Schenk's  CHIEF  method 
at  the  Helmholtz  eigenfrequencies.  However,  Patel's15  like  use  of  CHIEF  as 
part  of  his  analysis  of  scattering  from  elastic  structures  is  unnecessary 
by  the  above  argument  since  he  ignores  the  fact  that  the  in-vaauo  structural 
formulation  represents  interior-problem  enough  for  the  external  fluid. 


11 
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The  modal  fluid  model  developed  here  is  essentially  an  extension 
of  those  by  Sandman  to  include  the  effect  of  nonzero  cap  motion.  Thus, 
inversion  of  the  fluid  system  derived  in  Chapter  II  would  yield  modal 
impedances  immediately  suitable  for  comparison  to  those  of  Ref.  10. 
Although  such  results  were  generated  as  part  of  the  benchmark  shake- 
down  of  the  present  fluid  model,  for  the  sake  of  conciseness  we  do  not 
report  them  here. 


B.l.b  Asymptotic  Models 

Earlier  we  stated  that  the  total  cost  for 

solving  the  two  coupled  systems  for  a  given  arbitrary  frequency  is 

sometimes  dominated  by  the  generation  of  the  fluid  system.  Add  to  this 

the  need  to  generate  a  new  fluid  system  for  each  new  frequency,  and  the 

incentive  for  a  rational  substitute  for  the  Helmholtz  equation  formulation 

16 

becomes  clear.  Geers'  Doubly-Asymptotic  Approximation  (DAA)  requires 
in  effect  that  the  fluid  system  be  generated  for  only  one  frequency:  the 
zero  frequency  which  yields  the  apparent  mass  matrix  for  the  given  geometry. 
The  method  has  been  refined  to  give  the  exact  fluid  impedance  not  only  at 

the  low  and  high  frequency  limits,  but  also  at  an  additional  midrange  value 

17  18  1  19 

of  choice.  DAA  has  been  applied  in  a  number  of  studies  '  '  which 

gauge  its  performance  by  comparing  to  exact,  analytically  available  results. 

20 

In  a  more  general  application,  Rantlet  and  DiMaggio  have  used  DAA  to  model 
radiation  loading  for  transient  response  to  shock  of  a  finite  cylindrical 
shell  for  which  the  alternative  would  be  the  more  expensive  Helmholtz  integral 
formulation. 

In  Ref.  1  DAA  was  used  as  part  of  an  iterative  technique  to  compute 
the  exact  structural  response  corresponding  to  an  exact,  Helmholtz-integral 
fluid  model.  DAA's  purpose  was  to  stabilize  iterations  whose  convergence 
was  otherwise  either  slow  or  impossible  given  the  schemes  devised.  Con¬ 
vergence  for  the  techniques  here  met  with  no  such  difficulties  and  a  similar 
recourse  to  DAA  was  therefore  not  necessary. 


4 
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B. 2  Interation  Models  for  Finite  Cylinders  using  "Exact" 

Two-  and  Three-Dimensional  Fluid  Formulations 

21 

Chen  and  Schweikert  have  applied  a  general  finite- 

element  methodology  to  a  cylindrical  shell  and  have  modeled  the  three- 

dimensional  fluid  interface  as  a  superposition  of  sources  of  strengths 

unknown  a  priori.  An  internal  stiffening  structure  was  also  considered. 

In  accordance  with  the  earlier  discussion  here,  they  met  with  no  matrix 

22 

inversion  problems.  Baron,  et.  al.  used  a  finite-difference  description 
of  a  solid  elastic  cylinder  by  dividing  it  into  bands  and  assuming  a  modal 
expansion  in  the  circumferential  direction.  The  two  cap  surfaces  were 
similarly  divided  into  concentric  bands.  The  three-dimensional  fluid  was 
modeled  using  the  Helmholtz  integral.  The  final  coupled  system  was  obtained 
by  equating  the  surface-normal  stress  and  displacement  fields  to  the  ex¬ 
ternal  fluid  pressure  and  particle  displacement,  respectively.  A  later 

model  which  includes  the  possibility  of  an  internal  void  is  in  fact  a 

23 

fluid- loaded  thick-shell  theory  for  the  finite  cylinder.  Later  yet,  an 

orthotropic  shell  model  was  developed  to  simulate  stiffening-frame  effects, 

24 

and  elastic  end  caps  were  incorporated. 

25 

Crighton's  analysis  of  acoustic  scattering  by  a  finite  cylindrical 
bar  not  only  models  the  fluid  as  two-dimensional,  but  has  also  adopted  a 
strip- theory ,  or  locally-reacting,  radiation  approximation  to  eliminate 
the  cumbersome  fluid  loading  integral.  Such  simplifications  allowed  him 
to  arrive  at  a  simple  and  useful  theoretical  result  for  the  ratio  of  res¬ 
onant  to  nonresonant  far-field  levels. 


rroughs  ,  and  Cole  have  also  used  two-dimensional  fluid  models 

in  th _ _  lalyses  of  force-driven  framed  cylindrical  shells.  Both  studies 

were  in  go  agreement  with  experimental  data.  Cole  models  each  stiffening 

frame  as  a  Lng  of  finite  thickness  and  width,  and  calculates  in-and  out- 

„  .  .  28 

of-plane  displacements  due  to  shell-reaction  loads  by  means  of  Harari's 
equations . 


9 

5 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC,  • 


The  familiar  far-field  approximation  for  the  radical  in  the  exponent 
is  then  applied. 


2  —  2  ( 
/l+r  +2rcosa+ (z-z)  =  R <  1  + 


( ,  rcosa-zz  ) 

1  +  2  / 

(II.A-23) 

R  ’ 

which  may  be  further  simplified  to  just  R  for  the  radical  in  the  denominator 
of  Eq.  II.A-22,  which  becomes 


ikaR  ikaR  ...  2a  -  .  mira  -  -ika(z/R)z 

e  .  (1)  .  .  e  ,  , . n+1  /  dz  sin  - — (z+L/2a)e 

— r —  f  (r,z)  =  — - —  -  (-1)  /  L 

R  nm  R  / 


-L/2a 


’  J  dotcosna  e 
0 


ika (r/R) cosa 


(II.A-24) 


The  z  integral  is  easily  evaluated.  For  the  a  integral  one  recognizes  that 


rn  r  if, 

J  ( t)  =  -  I  dacosna  e 

n  tt  J 


so  that  finally 


f  (0)  =  -  — — r -  J  ( kasindi ) 

nm  2  n 


J  imiT/2  sin  [ ( imra/L-kacosiji }b/2a 
I  f  mira/L-kacosifj} 


- i mii/2  sin  f  {mna/I'.+kacosi^)L/2a] 


{ mu  a/L+kacosiL } 


(II.A-25) 


Applying  a  similar  reasoning,  one  finds  that 


t  ^  (,(,)  =  kasinij,  /  J  (kasi 

nm  1  4  )  n+1 


iniji)  +J ,  .  (kasiniji)| 


{i  mit  /  2  sin  [  { mna/L-kacosi[i)L/2a  | 
[miia/h-kacosi|j) 


-imti/2  sin  [  {nnra/L+kacos'M  L/2a]  | 
{mna/h+kacosijd  f 


( T I .A- 26) 
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ika  pc  T. ,  f  ^  V 

m=l  nm  nm 


+  F  f(2)  p 
m=l  nm  nm 


+  ika  pc  f(3)  V  L 
n  n 


—  T.  „  cosnT<  ,  C 

R  n=0  J  c(4)  L 

+  f  p 
n  n 


(II.A-21) 


(5)  CR 

4-  L>;a  pc  f  V 
n  n 


(6)  °R 
+  f  p 
n  n 


where  we  have  analytically  extracted  out  the  free-field  asymptotic  factor  e 

(R  nondimensional ) of  the  three  dimensional  pressure  field.  The  factor  multi- 
lkaR 

plying  e  /R  on  the  right-hand  side  of  Eq.  II.A-21  is  the  directivity 

function  for  the  vibrating  shell  with  end  caps  as  a  complex  radiator  (the 

angles  <£ ,  ip  are  defined  in  Fig.  II. 1).  Thus,  the  quantities  f^l...,f^ 

nm  n 

are  functions  of  (r/R,z/R)  (sin  ip ,  cos  ij>)  available  analytically  by  approxi¬ 
mating  the  integrals  in  Eq.  II.A-7  resulting  from  substitution  of  Eqs.II.A-9a 
through  II.A-9f.  Each  of  their  contribution  to  the  total  field  has  a  simple 

physical  interpretation.  For  example,  f^^,  f^2  represent  modal  monopole 

nm  nm 

and  dipole  contributions  over  the  shell  surface  since  respectively  they 

(3) 

appear  multiplying  V  and  p  .  Similarly,  f  represents  the  directivity 

nm  nm  ,  . .  n 

(4) 

contribution  of  a  left  cap  monopole,  and  f  that  of  a  right  cap  dipole, 

n 

etc.  The  expressions  given  below  clearly  show  all  these  features.  To  cal¬ 


culate  f 


(ip)  ,  one  starts  with  its  exact,  near-field,  form: 


ikaR  . , . 

e _  c  (1)  , 

— r—  f  (r 

R  nm 


,  3  C  j-  irma,-  . 

,z)  =  -  —-  /  dz  sin - (z+I,/2a) 

4  it  J  L 


-b/2a 


(II.A-22) 


'  2 (-  1 ) 


n  /  dotcosnra  e 


/  2  -  2 
ika/l+r  i2rcosa+ ( z-z) 


2 


2 
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and 


L/a  .  imraz  f  cos  \  (m-m)  7T  (1-az/L)  ] -cos  (m-m)  tt/2 
»  sin  \  z 

2  it  LI  m-m 


cos  [  (m+m)  7t  (1-az/L)  1-cos  (m+m)  tt/2 
m+m 


L/a  nnraz  j  sin  [  (m-m)  tt  { 1-az/L)  ] -sin  (ra-m)  tt/2 

+  -  cos  -  \  - - - 

2  it  L  l  m-m 


sin  [  (m+m)  tt  (1-az/L)  ]  -  sin  (m+m)  tt/2  J  for  m  _/ 
m+m  /  r 


m 


(II. A- 


where  again,  here 


al  (z)  = 
In 


2  (-1) 


A 2+z 2  0 


TT 

—  J~  dacosnaexp{ika>^+z2  X+gcosg 


(II. A- 


/  1+  fcosa 


a  !2n(Z) 


=  2(-l)n  f 

~  {2+z2}3/2  J 


dacosna  (1+cosa)  [1-ika 


A+z2  A 


l+^cosoJ 


{ 1+Ccosa} 


3/2 


with 

E  =  2/(2+z2) 


■  expfika 


v4+z2  A 


l+£cosa} 


(II. A- 


2.  Calculation  of  the  Acoustic  Field 


pressure  p 


ff 


Eq.  II.A-7  is  now  used  to  calculate  the  acoustic  far-field 
defined  as  (with  R=/ r^+z^ ) 


P 


ff 


lim 

/  2 


2 

+  7  -  vrr> 


n^O 


cosrnip  (r,z) 
n 


lim 

R-X» 


I.  cosntjip  (r,z) 
n=0  n 


18b) 


19a) 


19b) 


(II.A-20) 
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muaz  i  ,  .  ,,  2az,  .  .  2mir, 

— — —  <  cos  (mit/2)  [mir(l  -  — — )  +  sm  — — 

■  /  .  r  2muaz  .  ,,m,  I 

+  sin  (imr/2)  [cos  — - (-1)  ] 

muaz  (  ,  2muaz  ,  ,  ,m, 

— — —  <cos(mir/2)  [cos  — — —  -  (-1)  ] 


+sin (mir/2)  [mn(l  -  ~£Z)  -  sin  ?.m^.az ]  j  for  m  =  m 


muaz 


/  ,mu,  rsi 
Scos(— )  [— 


( II . A-17a) 


sin  (m-in)  ir/2-sin  (m-in)  itza/L 
m-m 


_  sin  (m+m)  ir/2-sin  (m+m)  uza/L 
m+m 


.  .m-rr.  , cos  (m-in)  ir/2-cos  (m-m)  nza/L 

+sm  (— )  [ - — - 

2  m-m 


.cos (m+m) n/2-cos (m+m) uza/L.  ) 
” '  m+m  J  | 

mirza  (  ,  .cos (m-m) u/2-cos (m-m) uza/L 

n  —  |  cos  (mu/2)  [ - — - 

,  cos  (m+m)  u/2-cos  (m+m)  uza/L. 


m+m 


.sin (m-m) u/2-sin (m-m) uza/L 

-sin (mu/2  - ^ 


+  sin (m+m) u/2-sin (m+m) uza/L,)  fpr  m  ^  ^ 


m+m 


(II.A-17b) 


{< 


muaz  I  2miraz  .  ..ml 

in  — : —  (cos  — - (-1)  f 


os 


muaz  (  2az,  .  2muaz\ 

— - —  <  (1 - — )  mu  +  sm  — = - >  for 


m  =  m  (II.A-18a) 
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the  degree  of  computation  is  further  reduced.  Finally,  all  entries  for 

which  m+m  is  not  an  even  integer  may  be  shown  to  be  identically  zero.  It 

2 

turns  out  therefore,  that  instead  of  M  triple  integrals  for  the  MxM  sub¬ 
matrix,  we  actually  have  only  M/2* (M/2+1)  double  integrals  for  M  even,  or 
2 

(M+l)  /4  double  integrals  for  M  odd.  Thus,  for  example,  if  M=30  the 

number  of  double  integrals  needed  for  the  MxM  submatrix  becomes  15 (15+1) =240 

(not  900).  If  M=31,  then  (31+1)2/4  =  256  double  integrals  are  required. 

2 

For  large  M  therefore,  roughly  M  /4  double  integrals  are  needed  to  generate 

most  of  the  fluid  elements  by  the  present  modal  method.  Since  finite  dif- 

2  2 

ferences  would  require  M  (or  perhaps  M  /2  due  to  reciprocity) ,  the  present 
method  represents  a  net  savings  for  this  part  of  the  calculation.  However, 
the  matrix  entries  for  the  bottom  two  rows  and  rightmost  two  columns  of  [C] 
and  [D] ,  require  the  evaluation  of  roughly  3M  triple  integrals  for  n=0  or 
2M  for  n  >  1  by  the  modal  method  and  the  same  number  of  double  integrals  by 
finite  differences,  so  that,  calling  N  the  number  of  mesh  points  in  the  out¬ 
most  of  our  triple  integrals,  it  follows  that  both  techniques  are  roughly 

2 

equivalent  in  cost  when  3MN  =  3/4M  . 

The  double-integral  versions  of  D^1-  turn  out  to  be: 

nmm  nmm 


ca! 

nmm 


/  dz  {(T  +T  )  [a2I  (z)]+T0[a2I„  (z+L/2a)]} 
J  1  3  2n  z  zn 


(II.A-15a) 


[1_(_1)ITH'm]  f  dz  {(T  +T  )  [al  (z)]+T  [al  (z+L/2a)]} 
nmm  2  it  J  1  3  In  2  In 

0  , 


( II . A-15b) 


where  the  quantities  T^ ,  T2,  T3  are  functions  of  m,m  and  z: 


L/a,  ,  ,.m,  ,  rmraz  L/a  rmraz  c 

T  =  [l-(-l  sin  —  +  ~  cos  —7—  for  m  =  m 

1  4  Trm  L  4  L 


(II.A-16a) 


L/a  .  muaz  (  (-1)  -cos  (m-m)-rr/2  (-1)  ^ -cos  (m+in)n/2 

=  -  sin  -  <  - :: - +  - r= - — 


rn  =  -  q  I  n  — -  < - - - - -  + - - - - - XT 

1  2TT  L  )  m-m  m+m 


L/a  mnaz  )  sin  (m-m)Tr/2  sin  (m+m)iT/2  ( 
-  ■  -  r  cos  — —  \ = -  -  - r= -  / 


for  m  ¥  m;  (II.A-16b) 
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CL  CR 

For  n  >  1  the  present  model  has  V  =  V  =  0;  so  that  the  right  side  of 
—  n  n 

Eq.  II.A-13  becomes  then 


(1)  „(D 


nil 


3(1) 

n21 


3d) 

nMl 


nl2 


3(2)  -D(2) 

nm=l  nm=2 


3(1) 

nlM 


3(1) 

n2M 


3(1) 

nMM 


(-1)M+1d(2)m 

nm=M 


nm=l 


nm=2 


V 


nm=M' 


(II.A-14) 


where  the  last  two  columns  have  been  deleted  from  the  [D]  matrix  multiplying 
{v},  which  has  been  shortened  accordingly.  If  {V}  short  above  were  known  so 
that  the  above  product  could  be  performed,  thus  yielding  a  M+2-element  vector, 
the  system  governing  the  {p}  vector  would  still  have  rank  M+2;  its  solution 

would  give  the  shell  modal  fluid  loading  p  p  w  plus  the  fluid  cap 

p  ~  nm-1  nm=M 

W  L-p  , 

pressures  p^  ,  Pn  •  From  the  acoustic  viewpoint,  these  cap  pressures,  whether 
or  not  there  is  cap  motion,  represent  the  modal  strength  of  radiating  cap  di¬ 
poles. 

All  constants  in  both  [C]  and  [D]  matrices  require  the  evaluation  of 
a  triple  integral.  Thus,  the  present  modal  formulation  initially  does  not 
appear  economically  advantageous  compared  to  a  similar  finite-element  or 
finite  difference  calculation,  e.g.,  that  in  Ref.  22  where  only  double  in¬ 
tegrals  are  required.  However,  we  now  show  that  the  bulk  of  the  matrix 
elements,  namely  those  composing  the  upper  left  MxM  submatrix  in  both  [C] 
and  [D ]  matrices  which  represent  the  coupling  of  m,m  modes  on  the  cylindrical 
part  of  the  surface,  may  be  recast  in  terms  of  double  integrals.  Moreover, 
these  submatrices  are  symmetric,  e.g.  in  Eq.  II.A-11  C^~  =  ,  so  that 
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Finally,  if  the  modal  sum  is  everywhere  truncated  to  ,  M  being 
an  integer  of  sufficient  size,  and  advantage  is  taken  of  certain  symmetry 

relations  between  some  of  the  constants  (e.g.  ,  it  turns  out  that  D^  =  (-l)m+1 

( 2 )  ( 3 }  m ] i  ( 2 )  nm 

D  ;  C  =(-l)  .  C  ) ,  the  resulting  set  of  fluid  equations  may  then  be 

nm  nm  nm 

cast  in  the  following  matrix  form: 


L/2a+C(^ 

nil 

C(1) 

nl2’  ‘ ' 

c(1) 

nlM 

!  c<2\ 
i  nm=l 

i 

C(2,1 

C(1) 

n21 

L/2a+Cn22“' 

c(1) 

n2M 

i  r(2) 

!  nm=2 

i 

-c(2) 

nm=2 

C  ^ 
nMl 

... 

L/2“*cnm 

i 

i 

i 

!  c(2) 
i  nm=M 

i 

<-dm+1c(2)m 

nm=M 

.. 

£(2\ 

nm=l 

c<2>  ... 

nm=2 

c  (2) 

nm=M 

"T - 

:  1/2 

i 

i 

£(3) 

n 

.  ...  . 

_c{2> 

nm=2 

nm=M 

I 

!  -a(3) 

i  n 

■ 

1/2 

D(1>  ... 

nl2 

D(1) 

nlM 

i  o(2\ 

|  nm=l 

i 

°(2\ 

nm=l 

D(1) 

n22 

•  • 

D(1) 

n2M 

!  d(2), 
i  nm=2 

i 

f 

-d(2) 

nm=2 

•  • 

D(1> 

nMM 

i 

|  d(2) 

,  nm=M 

f 

- 

<-i>KtV2’ 

nm=M 

_ _ _ 

i 


i(2)1 

nm=l 

d(2) 

nm=2 

d(2) 

nm=M 

I 

1 

1 

-f— 

d(3) 

n 

fi(4) 

n 

) (2) 

-d(2) 

(-DM+1d(2)m 

r 

5(4)  i 

d(3) 

nm=l 

nm=2 

nm=M 

! 

n  i 

1 

n 

^nm=l 


p 

nm=2 


^nm=M 


V 


nm=l 


V 


nm=2 


V 


nm=M 


Vn=0 


n=0  J 


(II.A-13) 
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L/2a 


„(1)  ika  f  , 

D  ~  —  /  dzsi 

nmm  2tt  j 


-L/2a 
L/2a 


L/2a 

/ 

-L/2a 


sin  52®  (z+L/2a)  i  d£sin  ~  (z+L/2a)  [alln (q, r=l)  ] 


,(1) 


nmm 


/ 


-L/2a 
L/2a 


L/2a 

f 6 

-L/2a 


dzsin^^-  (z+L/2a)  -f  dzsin^^-  (z+L/2a)  [a2I  ( £  ;r=l)] 
L  /  L  2n  1 


,  ( 2) 


nm 


ika 

2tt 


/ 


dzsin  — ^ —  (z+L/2a)  J  drr[al  (£  ;r=l)] 


1 

r  - 


L  - »  - 3n  2 

-L/2a  0 


L/ 2a 


.(2) 


nm 


/ 


dzsin  (z+L/2a)  •  (z+L/2a)  J  Arr  [a3I4n  (?2;  r=l)  ] 


-L/2a 

L/2a 


0 


(3) 


nm 


ika 

2tt 


/ 


dzsin 


mira 


-L/2a 
L/  2a 


/• 


L  (z+L/2a)  j  drr[al5n(?3;r=l)  ] 


C(3)  =  -  —  f  dzsin  ^  (z+L/2a)  •  (z-L/2a)  f  d??[a3I 
nm  2tt  L  J 


6n(?3;r=l)]  (II . A- 


-L/2a 


The  constants  CK  L  and  l  have  been  initially  defined  in  terms  of  Cauchy- 
nmm  nmm 


Principal  Value  integrals  which  omit  the  point  z=z  for  which  £^=1  and  for  which 
2, 


al  and  a  I  are  therefore  infinite.  (Actually,  it  may  be  shown  that  the 
In  2n 

resulting  singularity  is  of  the  log  type,  and  thus  integrable  so  that  the 
Cauchy-Principal  Value  designation  may  be  lifted.  Such  a  result  is  consistent 
with  the  discussion  in  Ref.  34.) 


This  procedure  may  be  repeated  for  field  points  on  the  left  and  right 


caps.  For  the  point  on  the  left  cap,  for  example,  the  left-hand  side  of 


Eq.  II.A-7  becomes  p  and  the  second  integral  on  the  right-hand  side  must  now 

n 


be  written  as 


f 


i.e.,  we  return  to  Eq.  II.A-7  and  let  z=-L/2a  ;  multiply  then  by  rdr  and 
integrate  both  sides  of  Eq.  II.A-7  from  0  to  1 .  Constants  similar  to  those  in 
Eqs .  II .A- 11  and  II.A-12  may  be  defined;  these  are  given  in  Appendix  II. A. (1).  A 
similar  calculation  is  performed  for  the  right  cap. 


12) 
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P  (r,-L/2a)=P  (-L/2a)  =  p  ,  a  constant  function  of  end-cap  (II.A-9c) 

n  radial  coordinate  r 


p  (r,Ij/2a)  =P  (L/2 a)  =  p  ,  a  constant  function  of  end-cap 
n  radial  coordinate  r 


(II . A-9d) 


V  (r,-L/2a)  = 
n 


V  (r,L/2a)  = 
n 


V  ,  a  constant  function  of  end  cap  radial 
coordinate  for  n=0 


0  for  n  >  1 


V  ,  a  constant  function  of  end  cap  radial 
coordinate  for  n=0 


0  for  n  >  1 


(II.A-9e) 


( II . A-9f ) 


Similar  expansions  may  be  used  for  the  field  pressure  pn(r,z)  on  the  left- 
hand  side  of  Eq.  II. A- 7  as  the  point  approaches  the  cylinder  surface;  e.g.  ,  for 
a  field  point  (r,z)  now  on  the  cylindrical  shell  surface  (l,z)  we  have 


oo  mTTa  00  I 

p  ~  sin  —■  (z+L/2a)  =  < ikapcV  - 

m=l  nm  L  m=l  l  nm 


'  ^  / 


—  .  mTra  —  _  ,  „  . 

dz  srn  — —  (z+L/2a) 
Li 


-L/2a 


[aIm(?i ?r=1)1  +  p 


(L/2a 

/ 

-L/2c 


dz  sin—™  (z+b/2a) [a2I2n (Cj ; r=l) ] 


(II.A-10) 


where  now  £  =  2/[2+(z-z)  ],  and  are  similarly  evaluated  for  r=l . 

Multiplying  the  above  by  sinmnf./L (z+L/2a)  and  integrating  from  -L/2a  to 
+L/2a  with  respect  to  z  yields 

C 

oo  (1)  oo  (1)  (2)  L 

L/2a  p  =  pc  Dv  1  V  ~  -  Cv  L  p  .  +  pcD'  V 

nm  m=l  nmm  nm  m=l  nmm  nm  nm  n 

(II. 

-  c  (2»  P°L  ♦  pcd<3»  -  c<3>  pS 

nm  n  nm  n  nm  n 


(II.A-11) 


where 


11 
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Pn(r'z)  =  ~  4u 


L/2a 

/ 

-L/2a 


dz{ ikapcV^ (1 ,z) [al(r,l,z,z)l 


+p  (l,z)  [a  I  (r,l,z,z)]} 
n  2n 


e 

4tt 


j. 

f 


drr{ikapcV^(r,-L/2a) [al (r , r , z , -L/2a) ] 


(II.A-7) 


+(z+L/2a)p  (r,-L/2a) [a  I  (r , r ,z, -L/2a) ] 
n  4n 


£ 

4tt 


j. 

/ 


drr{ikapcV  (r,L/2a) [al_  (r,r,z,L/2a) ] 
n  5n 


where, for  example, 


- (z-L/2a) p  (r ,L/2a) [a  I  (r , r ,z ,L/2a) ] } 
n  bn 


alln(r,l,z,z)  = 


2  (- 


t  2  2 

ir  ika/l+r  +(z-z)  /l+f^cosa 

dacosnae 


/l+r2+ (z-z) 2  0 


/  1+E^cosa 


(II.A-8) 


with 


= 


2r 


„  2  ,  -.  2 
1+r  +(z-z) 


2  3 

The  rest  of  the  expressions  for  a  I  ,  a  1^  are  given  in  Appendix  II. A. 

In  Eqs.  II.A-7  and  II.A-8,  all  spatial  variables  have  been  normalized  by  the 

cylinder  radius  a.  V^(l,z)  and  pn(l,z)  stand  for  the  normal  velocity  and  pressure 

on  the  cylindrical  part  of  the  boundary;  V  (r,-L/2a)  and  p  (r,-L/2a)  stand  for 

n  n 

corresponding  quantities  at  the  left  cap;  and  V^(r,L/2a)  and  p^(r,L/2a)  stand  for 
those  at  the  right  cap.  We  now  assume  the  following  expansions: 


V  (1,5)  =  E ,  V  sin  5—  (z+L/2a) 
n  m=l  nm  L 


(II . A-9a) 


p  (l,z)  =  T  p  sin  (z+L/2a) 

n  m=l  nm  L 


(II .A- 9b) 
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A.  The  Fluid  System 


1 .  Radiation  Loading  System 


In  a  fluid  the  Helmholtz  integral  relates  the  velocity  potential 4> 
at  a  field  point  to  boundary  values  of  potential  and  velocity,  both  initially 
unknown  in  an  interaction  problem: 


-IS 


~  34>  3G  . 

dS  /  G - 4> 

I  3  n  3  ri 


(II.A-1) 


where  e  =  1,2  for  points  in  the  fluid  and  on  the  boundary,  respectively.  G  is 
the  free-field  Green's  function, 


(II.A-2) 


where  k  is  the  acoustic  wavenumber  w/c,  and  R,  the  distance  between  field  point 
(r,(|>,z)  and  boundary  point  (r,a,z),  is  given  by 


/  2  -2  -  2 

R  =  /r  +r  -2rrcos +  (z-z) 


(II.A-3) 


so  that 


for  points  on  the  right  cap  (+)  and  left  (-) . 


(II.A-4) 


I  —  for  points  on  the  curved,  longitudinal  part  of  the 
l  3r  cylindrical  surface. 

The  r,(f>,z  coordinate  system  is  shown  in  Fig.  II. 1;  (r,a,z)  correspond  to 

running  values  of  these  variables  in  the  integration. 

The  fluid  pressure  anywhere  may  be  found  from  <}>  from  the  unsteady  linearized 
Bernoulli  equation 


p  =  “  p  aft  =  lwp$ 


(II.A-5) 


If  one  writes  p  in  terms  of  circumferential  n-modes  as 

00  (1  )  00  (  2) 
p(r,$,z)  =  p^  (r,z)cosn<|>  +  £  p^  (r,z)sinn<j> 

it  may  be  shown  that  both  and  p^^  satisfy  the  following  form  of  the 

Helmholtz  integral: 


(II.A-6) 


l 
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II.  FLUID  AND  STRUCTURAL  MODELS  FOR  A  FINITE,  FRAME-STIFFENED  CYLINDRICAL 

SHELL  KITH  END  CAPS 

As  earlier  pointed  out,  the  finite  cylindrical  geometry  is  nonseparable. 
Because  of  this,  the  modal  form  chosen  to  describe  pressure  and  displacement 
variations  along  the  axial,  or  curved,  part  of  the  surface  will  inevitably  be 
nonorthogonal  to  that  chosen  to  describe  these  same  variables  over  the  cap 
surfaces.  The  result  is  an  elaborate  set  of  linear  equations  which  for  both 
structure  and  fluid  couples  one  set  of  modes  to  the  other. 

Here  we  assume  that  the  simply-supported  set  of  sinmirz/L  functions 

adequately  describes  fluid  and  structural  quantities  over  the  longitudinal 

part  of  the  geometry — an  assumption  of  course  not  met  by  displacements  near 

the  ends  of  a  freely  suspended,  or  neutrally  bouyant,  cylinder  but  which 

nevertheless  has  proven  successful  in  predicting  acoustic  radiation  from 
27 

such  structures. 

Also,  because  the  analysis  is  here  ultimately  applied  to  a  case  for  which 
the  cylinder  radius  is  much  smaller  than  the  shortest  acoustic  wavelength 
considered,  only  the  radially  uniform  mode  is  used  for  the  caps,-  i.e.,  cap 
displacements  and  pressures  are  given  by  the  product  of  the  "unit"  mode  and 
the  constants  which  denote  their  values.  Should  future  need  dictate  a  more 
complete  cap  description,  an  expansion  such  as  that  of  Ref.  20,  Eq.  62  may  be 
used. 

In  Section  II. A. 1  the  Helmholtz  integral  is  used  to  develop  a  radiation 
loading  model  for  the  finite-cylindrical  geometry.  In  Section  II. A. 2  the 
fluid  system  so  generated  is  approximated  analytically  for  an  arbitrary  far- 
field  point.  The  radiated  pressure  thus  appears  as  a  function  of  surface 
pressure  and  displacement,  both  unknown  at  this  stage  in  the  analysis. 

Section  II. B  contains  the  analysis  for  the  dynamics  of  a  submerged  frame- 
stiffened  shell  with  end  caps.  The  otherwise  arbitrary  external  drive  is 
assumed  axisymmetric  in  the  circumferential  direction.  Together  with  the  fluid 
model  of  Section  II. A,  the  structural  system  developed  here  defines  the  fluid- 
structure  interaction  problem. 
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Computational  Fluid  Dynamics  almost  by  definition  deals  with  very 

large  linear  systems,  and  workers  in  that  field  have  often  resorted  to 

iterative  techniques  as  the  only  practical  means  of  solution.  An  example 

of  a  predictor-corrector  type  of  CFD  application  of  iteration  arises  in 

the  time-domain  description  of  the  fluid  medium  bounded  by  a  complex 

structure  in  impulsive  motion,  a  calculation  which  normally  requires  the 

development  of  the  solution  over  a  large-volume  grid  of  mesh  points  at 

each  time  step,  and  the  associated  enforcement  of  continuity  of  particle 
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velocity  at  each  instant  over  the  contact  surface.  Chang  and  Wang  have 
used  an  iterative  scheme  in  their  transient  boundary  condition  calculation. 

At  the  end  of  each  time  step  the  estimated  fluid  pressure  was  used  to  cal¬ 
culate  structural  displacements,  which  then  generally  did  not  agree  with 
those  for  the  fluid  particles  at  the  interface.  These  were  then  corrected 
and  the  fluid  medium  recalculated  to  yield  the  new  radiation  load  vector, 
etc. 

32 

Belytschko  has  discussed  the  role  of  iteration  in  the  implicit 

time  integration  of  system  equations  for  the  general  transient  fluid- 
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structure  interaction  problem.  Whitlow  and  Harris  have  done  essentially 
the  same  for  unsteady  transonic  flow  in  their  discussion  of  similar  numer¬ 
ical  techniques  commonly  used  in  that  field.  Transonic  flow  over  an  oscil¬ 
lating  airfoil  is  governed  by  a  nonlinear  potential  equation  which  gives  rise 
to  shocks  in  the  gas  medium;  since  the  positions  of  these  vary  with  time, 
a  fine  local  mesh  is  usually  required. 

As  will  be  discussed  in  Chapter  III,  in  developing  some  of  the 
iteration  techniques  here  we  have  emphasized  the  convergence  role  of 
system  subpartitioning  rather  than  the  application  of  methods  known  to 
accelerate  convergence. 
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The  present  structural  model  is  essentially  this  same  model  with 
rational  simplifications  to  reduce  the.  number  of  degrees  of  freedom 
matched  between  shell  and  frames  (here  only  radial  and  circumferential 
displacements  are  used  whereas  Cole  also  matches  axial  displacements  and 
out-of-plane  rotations) .  Fundamental  differences  between  the  model  here 
and  that  in  Ref.  27  are  therefore  really  confined  to  just  the  fluid- loading 
system,  for  which  the  present  study  may  be  viewed  as  providing  the  three- 
dimensional  extension. 

B. 3  Fluid- Structure  Interation  Studies  using  Iteration  Models 

g 

Chertock  has  used  an  iterative  technique  to  effectively 
invert  the  fluid  system  for  a  general  body  of  revolution,  in  essence,  to 
find  its  impedance  matrix.  As  discussed  earlier,  he  generally  found  that 
while  the  method  converges  easily  at  low  frequencies,  it  later  becomes  un¬ 
stable  with  increasing  frequency.  He  circumvented  this  difficulty  by  alter¬ 
ing  his  basic  scheme  to  use  an  average  value  of  previous  iterates  where 
before,  for  the  lower  frequencies,  he  had  used  only  the  most  recently  avail¬ 
able  value.  Ref.  29,  Chapter  ll’s  general  discussion  on  such  summability 
techniques  provides  the  rationale  behind  their  potential  stabilizing  effect. 
We  too  applied  this  basic  idea  in  our  iterative  treatment  of  the  fluid- 
loaded  driven  sphere.'*' 

For  the  low-frequency  limit,  where  the  externally  unbounded  fluid 
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medium  has  a  purely  inertial  effect  on  the  vibrating  surface,  Au-Yang 
has  used  an  iteration  algorithm  to  find  the  fluid-loaded  resonances  of  a 
cylindrical  shell.  He  starts  it  out  by  assuming  each  resonance's  value 
and  then  uses  the  fluid  system  to  calculate  the  fluid  loading  at  that 
frequency.  Next,  this  pressure  solution  is  introduced  into  the  otherwise 
homogeneous  structural  system,  which  yields  a  new  root,  and  so  on.  He 
found  the  technique  to  converge  nicely  for  most  of  the  frequencies  cal¬ 
culated. 


t 
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As  previously  discussed,  f  ( )  stands  for  the  contribution  of  dipoles 

nm 

on  the  cylindrical  part  of  the  boundary;  for  every  n-m  combination  we  note 

the  factor  sin ip  in  the  above  expression.  The  rest  of  the  modal  directivity 

contributions  f ^ ^  , . . . , f ^  turn  out  to  be 
n  n 


f  (3>  (ip)  = 
n 


ika  (L/2a)  cosiJj 


j. 

/dr  r  J 

n 


(kasintjjr) 


(II.A-27) 


(4)  ,,,  (-i)n+1  ,  .  ika  (L/2a)  cosijj  / 

(tf<)  =  -  - - -  kacosthe  T  J  < 


f  (tjj)  =  -  - - -  kacosthe 

n  2 


/  dr  r  J  (kasinipr) 

0  n  (II.A-28) 


.(5)  (-i)  -ika  (L/2a)  cosip 

fn  w  =  -  ~ e 


f*-r 


(kasintjjr) 


(II.A-29) 


,  (-i)n+1  ,  ,  -ika(b/2a)  cosiJj  f 

(\p)  =  - - -  kacosijie  r  / 

2  0 


dr  r  J  (kasint/jr) 
n 


(II. A- 30) 


Again, f^3^,  f^  account  for  monopoles  at  the  left  and  right  caps  (acting 
n  n  ( 4 )  (*6 ) 

only  for  n=0  in  the  present  model) ;  and  f^  and  f^  represent  dipoles  also 

at  the  caps  (hence  the  extra  cosil)  factor  for  f^,f^). 

n  n 

The  grouping  mira/L  +_  ka  cosijr  appearing  as  argument  of  sin{}/{)  in 

Eqs .  II.A-25,  26  to  a  large  extent  determines  the  far-field  effect  of  sources 

distributed  along  the  shell  part  of  the  cylinder.  Physically,  supersonic 

axial  modes  may  be  expected  to  radiate  efficiently  and  in  the  ray  direction  rjr* 

given  by  cosip*  =  1/M^,  where  denotes  axial  Mach  number  for  the  m  mode. 

Since  the  modal  solution  is  of  the  form  sin  miTaz/L  exp  (-iwt)  -exp  [-iw  (t-za/c^)  ]  , 

with  c  =wL/m7T,  it  follows  that  modal  Mach  numbers  M  are  given  by  wL/nrr  •  1  /c= 
m  m 

ka/ (mira/L)  so  that 


— r —  +  kacosib  =  ka  {77+  cosi|r} 
L  —  M  — 

m 


(II.A-31) 
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Since  the  sin  {}/{}  group  has  its  maximum  value  of  unity  at  the  null 
of  Eq.  II.A-31,  the  analytical  solution  embodies  the  physical  requirement 
of  preferred  radiation  in  the  theoretical  ray  direction.  Subsonic  axial 
modes  (for  which  Eq.  II.A-31  never  vanishes)  lack  a  clear  preferred  radia¬ 
tion  direction  and  in  fact  would  have  an  identically  zero  far-field  were 
the  cylinder  infinite  in  length.  Their  acoustic  effect  for  the  finite 
geometry  is  that  of  a  pair  of  sources,  one  for  each  end  of  the  cylinder, 
and  thus  a  contribution  which  may  be  expected  to  be  generally  low  in  comparison 
to  that  of  supersonic  sources  distributed  over  the  entire  surface,  all  reach¬ 
ing  the  far-field  since  their  Mach  lines  (acoustic  wavefronts)  never  intersect. 

Similar  statements  may  be  made  about  radiation  in  the  circumferential 
direction  and  the  mathematical  behavior  of  J^(kasiniJ;)  as  a  function  of  order 
and  argument.  It  may  be  easily  shown  that  the  rotary  speed  of  circumferential 
modes  satisfying  n>ka  is  respectively  supersonic  and  subsonic.  Subsonic  con¬ 
tributions  instead  of  being  attributed  to  finite-end  effects  are  due  now  to  the 

sources'  felt  acceleration  as  they  travel  in  a  circle.  Such  is  the  prediction 

36 

of  the  general  FFowcs-Williams  Hawkins  equation  which,  incidentally,  also 
provides  the  theoretical  foundation  for  Gutin  noise  —  the  sound  radiated  by 
steady  blade  forces  in  a  subsonic  propeller.  Acoustic  radiation  from  the  sub¬ 
sonic  n>ka  circumferential  modes  is  here  really  due  to  the  same  mechanism  of 
acceleration  due  to  rotary  motion. 

An  important  feature  of  the  preceding  far-field  calculation  is  that  it 
accounts  for  possible  source  noncompactness  in  both  circumferential  and  axial 
directions.  This  is  necessary  because  in  Chapter  IV  the  model  will  be  applied 
to  a  very  high-aspect  ratio  cylinder  for  which  the  condition  kL>>l  is  satisfied 
over  much  of  the  frequency  range  investigated.  In  fact,  although  the  far- 
field  condition  kR>>l  is  always  strictly  satisfied,  the  ratio  R/L  of  field 
point  distance  to  cylinder  length  turns  out  to  be  only  about  4  (here  for 
convenience  R  is  dimensional) . 

However,  we  argue  that  by  keeping  just  the  zz  term  (and  no  higher- 
order  contributions)  in  Eq.  II.A-23  the  theory  adequately  accounts  for  phase 
differences  along  the  cylinder  axis,  or  equivalently  in  the  time  domain,  for 
the  z-dependence  of  retarded  times  t  =  t  -  R/c  =  T(z)for  sources  distributed 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


along  the  cylindrical  axis.  Given  that  the  approximation's  neglected 

2 

term  inside  the  parenthesis  of  Eq.  II.A-23  is  known  to  be  l/2(z/R)  ,  one 
may  readily  estimate  its  possible  impact  on  the  phase  information  for 
parameters  of  interest  here.  First,  we  note  that  axial  noncompactness 

will  have  greatest  effect  for  field  points  along  the  cylinder  axis  0=0 

o  -  2  -  ,  - 

or  180  .  For  these,  zz/R  =cos0 . z/R=z/R,  so  that  the  percent  error  of 
the  neglected  term  is  1/2 (z/R).  Thus,  for  the  Chapter  IV  values  of 
Zmax=  k/2=9.58  ft  and  aR=45  ft' ,  the  upper  bound  of  the  phase  error  is  about 
12%,  which  should  not  seriously  affect  the  character  of  signal  cancellation 
or  reinforcement  at  these  two  far-field  points. 

As  explained  immediately  after  Eq.  II.A-23,  the  radical  which  appears 
in  all  denominators  has  been  even  more  simply  approximated  to  yield  the 
three-dimensional  spreading  factor  1/R.  For  the  cylinder  of  interest  here 
the  greatest  "true"  source-field  distance  is  /45^+ (9 . 58) 2  =  46,  which 
implies  that  in  using  R=45  ft  as  measured  from  the  cylinder's  centroid  an 
error  of  less  than  2.3%  (or  .2  dB)  is  incurred  in  so  estimating  the  level 
for  that  "worst"  source. 

B.  The  Structural  System 

In  Section  A  we  have  derived  the  fluid  system  modal  equations.  Here 
we  develop  similar  equations  for  the  structural  system,  i.e.,  for  a  framed, 
simply-supported,  finite  cylindrical  shell  with  end  caps.  The  shell  is 
submerged  and  thus  loaded  by  the  fluid  pressure  vector  of  Section  A,  which 
we  recall  provides  for  radiation  loading  at  the  cap  ends.  Consistent  with 
the  analysis  in  Section  A  we  assume  that  cap  motion  takes  plase  only  for 
n=0,  so  that  as  earlier  stated,  computed  cap  pressures  for  n  1  are  due 
only  to  diffraction  from  the  cylindrical  part  of  the  shell  surface  and  may 
be  expected  therefore  to  yield  a  far- field  on-axis  level  of  modal  radiated 
pressure  significantly  below  that  for  n=0.  The  system  is  driven  by  a  force 
distribution  applied  either  at  one  of  the  frames  or  directly  anywhere  on 
the  shell  skin.  The  structure  may  also  be  driven  axially  by  a  load  acting 
at  the  end  cap.  Shell  frames,  caps, and  the  applied  drives  are  all  assumed 
circumferentially  axisymmetric  and  thus  describable  by  the  cosn0  series. 
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Sections  B.1,2  develop  the  model  for  an  unframed  shell  with  end  caps. 
Section  B.3  covers  the  inclusion  of  stiffening  frames. 

Methods  for  solving  the  structural-fluid  coupled  systems  will  be 
discussed  in  Chapter  III. 

1.  Unframed,  Capped  Cylindrical  Shell  with  Fluid  Loading 

The  equations  of  motion  for  a  finite,  simply-supported  cylindrical 
shell  appear  in  a  number  of  standard  references,  e.g. ,  Ref.  35. 


L^u+I^q+L^w  =  - 


2  2 
a  (1-y) 

Eh 


i  C  r  C 

<  pz  (4>)6<z+L/2a)-pz  (<|>)6(z-L/2a)S  (II.B-la) 


L  u+L.q+L-W  =  0  (II.B-lb) 

4  5  6 

2  2 

V+L8q+L9W  =  "  "hT"  /  <i>)  “Pr  ( z  »  4>  ^  }  (II.B-lc) 

where  u,  q,  w  stand  respectively  for  axial,  tangential  (or  circumferential), 

and  radial  displacements;  the  shell  has  radius  a,  Poisson  ratio  v,  thickness  h, 

and  Young’s  modulus  E.  Figure  II. 2  shows  the  generic  geometry.  The  effect 

of  the  caps  on  the  shell  has  been  modeled  as  a  set  of  (unknown)  reaction 
C  C 

pressures  p  in  Eq.  II.B-la;  shortly,  p  L'R  will  be  expressed  in  terms 

z  z 

of  cap  displacements  through  the  dynamic  equation  for  cap  motion.  F^  and 
p^  are  respectively  the  externally  applied  drive  and  radiation  loading  on 
the  axial  surface  of  the  cylinder.  The  operators  L^,...,L^  given  in  Appen¬ 
dix  II. B  correspond  to  the  Donnell-Mushtari  thin-shell  theory. 

For  a  simply  supported  cylinder,  we  let 


u  =  £  £  u 
n  m  nm 

q  =  H  q 

n  m  nm 


cosnA  cos  ^7-^  (z+L/2a) 
Li 

.  .  mu  a  , 

sinn<j>  sm  — —  (z+L/2a) 
Li 


(II.B-2a) 

(II.B-2b) 
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w  =  E  £  w  cosncj)  sin  (z+L/2a) 

n  m  nm  L 


(IT . B-2c) 


and  also 


C  C 

p  L' R ( <}> )  5  (z+L/2a)  =  E  T.  p  L,Rcosn(j>  cos  (z+L/2a) 

z  T  —  n  m  rz  L 

nm 


(II.B-3a) 


P„/F  =  T  E(p  ,F  )  cosnt}>  sin  (z+L/2a) 

x  x  ii  ill  x  x 

nm  nm 


(II . B-3b) 


Substituting  Eqs.  II.B-2,3  into  II.B-1  one  obtains 


A3 


nm 


nm 


nm 


Cb  °L 
P  -  p 
Znm  Znm 


2  ..  2\ 
a  ( 1— \>  ) 

Eh 


(II.B-4) 


nm  nm 


where  the  quantities  A^  through  Ag  are  given  in  Appendix  II.B(l). 

2.  Equation  for  Cap  Motion 

Each  of  the  two  end  caps  is  assumed  to  have  mass  M  and  to  move 

c 

rigidly  in  the  axial  direction  under  cap  fluid  loading,  drive  and  shell 
reaction  forces.  The  fluid  loading  on  the  caps  is  assumed  evenly  distri¬ 
buted  and  thus  constant  over  the  radius  as  measured  from  the  caps'  center, 
whereas  shell  reaction  forces  are  taken  as  acting  only  at  the  shell-cap 
circumferential  interface.  For  example,  for  the  left  cap  at  z=-L/2a. 


L  2 
-itoM  V  =  a 
c 


1 

/ 


C  V  C  , 

dr  r<p  L  6(r-l)  +  — — -  +  p  ' 

*  z  r  z 

n=0  o  n=0 


(II.B-5) 


where  we  have  accounted  for  the  possibility  of  an  externally  applied  drive 

CT  f  cr 

at  radial  position  r  ;  p  b'r  denotes  fluid  loading,  and  p  t  the  reaction 


o  z 


n=0 


n=0 
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with  the  shell  structure;  =  Pst?ra  ,  t  =  cap  thickness,  and  pg  =  density 
of  metal.  Eq.  II.B-5  thus  becomes 


C  C  C 

~cap  ,  L  ,  _  L  1  L, f 

pc  2  V  ^=p  +F  +  —  p 

n=0  z  z  2  z 

n=0  n  n 


(II.B-6) 


where  for  convenience  of  later  analysis  we  expressed  ZCap,  the  cap  reactive 
impedance  normalized  by  c,  as 


2Cdp  =  -  ika  (— — •)  (t/a)  it  (II.B-7) 

P 

where  p  is  cap  material  density.  A  similar  expression  may  be  derived  for 
CR 

V^_Q,  the  axial  velocity  for  the  right  cap.  Again,  the  present  model  allows 
for  a  finite  left  side  in  Eq.  II.B-6  only  for  n= 0.  As  discussed  in  Section  II. A, 
more  complex  cap  modes  could  be  modeled  to  match  higher  n  shell  modes  but  for 
simplicity  and  tractability  this  is  not  pursued  here.  Letting  now  z  =  +L/2a 
in  Eq.  II.B-2a  we  observe  that 


m=0 


(n=0) m 


( II . B-8a) 


Vn=0  m=0  U(n=0),m( 


(II . B-8b) 


where  U  =  -  iwu  which  from  the  system  in  Eq.  II.B-4  has  solution 
nm  nm 


2  2, 
a  (1-v  ) 

U  =iw  - - - 

nm  Eh 


!C  C 

B1(PzL  “  C  }  +  VFr  "pr  ]l 

nm  nm  nm  nm  \ 


(II.B-9) 


Constants  B^,  B^  are  given  in  Appendix  II.B(2) 


Since 


,  m 


2e  (-1) 

(5(z-L/2a)  =  Z  — -  cos  z+L/2a) 


(II.B-lOa) 
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<5(z+L/2a)  =  I  cos  —  (z+L/2a) 

in— 0  L/  a  L 


(Il.B-lOb) 


where 


1/2  for  m  =  0 


1  for  m  >  1 


Thus,  it  follows  that 


CT  2e  C 
L  m  L 

3z  L/a  Pz 

nm  n 


(II.B-lla) 


CD  2e  (-1)  CD 
R  m  R 


(II.B-llb) 


where  '  are  the  same  shell-cap  reaction  forces  appearing  in  the  equation 

for  the  cap  dynamics,  e.g. ,  Eq.  II.B-6.  These  ma^  be  eliminated  by  substituting 

Eqs .  II.B-11  into  II.B-9.  Using  Eq.  II.B-6  for  V  L  and  the  corresponding  right- 

n=0 

cap  equation  yields: 


2  2  -cap 

i  +  ima  (1-y  )  2pcZ  , 

Eh  L/a  m 


]C 

V  L 
n=0 


•  2,i  2.  -cap  C 

i»«_JirvJ.l2S5Z  B IV  / 

Eh  I  L/a  m=0  1  m  n=0 


cr  C 

2  L  2  m  R 

+  —7—  (  E  B  e  )  F  -  -4-{  E  B  (-l)e  )  F  K 

L/a  m=0  1  m  z.  ..  L/a  m=0  1  m  z, 

(n=0)  (n=0) 


C 

m  ,  R 


+  B.e  )p  L  -  Zn  B(-l)\  )p 

L/a  m=0  1  m  n  L/a  m=0  1  m  n 


Jil  ®oF_.  +  I.  B  p  \ 

T'-l  ’  r  m=l  3  *r  f 

nm  nm  I 


(II.B-12a) 
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2 2.  _  ~cap 

ioia  (1-v  )  _  2pc2 

Eh  L/a  nv 


B  e  )1  V  R  = 
i=0  1  m  J  n=0 


2  2 
iu>a (1-v 


SlzyjL  i  _2pc_2^r.  ,  E  B  e  )  v°L  -  -L-  (  r  B 

Eh  \  L/a  m=0  B1  Vn=0  L/a  n>=0  B1  X)  m)Pz(n=0) 


c  c 

+  77 —  (  T  B  e  )F  R  -  ~~  (  r  B  (-l)me  )pnL 
L/a  m=0  1  m  z.  ..  L/a  m=0  1  m  "n 

(n=0) 

C 

+  ~  (  I  B,  (-l)me  )P  R  +  E(-1)Vf  -  E  (-l)Vp 

L/a  m=0  1  m  n  m=l  3  r  m=l  3  r 

nm  nm 

(II.B-12b) 

°R  f  Ct  f  Cp  C 

where  we  have  simplified  the  notation  by  calling  p  '  p  '  simply  p  ,  p  — the 

2n  zn  n  n 

fluid  pressures  on  the  right  and  left  caps  as  denoted  in  Section  II. A;  and 

similarly  replaced  p  by  p  ,  the  modal  pressures  for  the  shell  part. 

nm 

From  the  system  in  Eq.  II.B-4  one  may  also  solve  for  v elocities 

V  =-imw  on  the  cylindrical  part  of  the  surface.  Substituting  all  we 
nm  nm 

know  so  far  into  that  solution  yields  that 


•  2  2,  (  2e  r  /  C  C  \ 

iaia  (1-v  )  1  ni  -cap  /„  L  ,  ,  ,m,T  r  1 

* - iii -  I  Vi  LpcZ  IVo-'-11  Vo) 

-  1-1)"fCr  )-i  (p  L  -(-l)’p  R)1  t  B I  (F  -p  )( 

'  2  (n=0)  z  (n=0)  2  '  “  n  'J  9  V  i 


(n=0) m 


(II.B-13) 


In  the  above  e  =  1  always  since  m  >  1.  For  n  >  1  the  cap  velocities 
C  m  — 

VnL'R  are  identically  zero  and  the  cap  fluid  loading  does  not  affect  the 
motion  of  the  cylindrical  surface.  Thus  for  n  >  1,  Eq.  II.B-13  simplifies  to 
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V 


n>l  ,m 


itoa2  (1-y2) 
Eh 


bq  (b  ~ 

9  r  nm 

nm 


(II. B- 


We  now  collect  the  results  in  Eqs.  II.B-12a,  b  and  B-13  for  n=0, 

and  II.B-14  for  n>Ll.  into  structural  systems  for  n=0,  >_1  respectively. 

2  2 

We  recognize  that  iwa  (1-v  )/Eh  =  iw(a/h) (p  c  (1-in))  (see  Appendix  II.B(l)) 

s  p 

and  define  the  following  three  constants 
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Then  for  n=0  we  have  that 


pc 


1 

0 

0 

• 

• 

0 

1 

1 

1 

m=l 

K1B7 

1 

1 

1 

K1B7= 

0 

1 

0 

• 

1 

1 

1 

1 

1 

vT2 

i 

1 

1 

1 

1 

-vT1 

0 

0 

• 

• 

0 

• 

1 

1 

1 

1 

1 

• 

1 

1 

1 

1 

| 

• 

• 

• 

• 

• 

0 

1 

1 

1 

1 

• 

1 

1 

1 

1 

• 

• 

• 

0 

1 

0 

1 

1 

1 

1 

• 

1 

1 

1 

1 

• 

0 

• 

• 

0 

1 

1 

1 

1 

1 

1 

1 

m=M 

K1B7 

» 

1 

1 

1 

1 

1 

(-DM+1kiB-m 

“1 

0 

• 

• 

• 

0 

1 

1 

1 

1+K1K2 

1 

1 

1 

1 

"Vi 

1 - 

■~T  ~ 

0 

• 

• 

• 

0 

1 

1 

-KK 

1 

1 

1  +  K  K 

V 


nm=l 


nm=2 


>m 


nm=M 

n=0 

-qr— 

tn=0 


14) 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC 


For  n  >_  1  the  structural  system  is  purely  diagonal 

n 

1  o 

l 

pc 

0  ‘1 


ifi (a/h) (p/p  ) (c/c  ) 
_ s  P 

/l-irt 


m=l 

b9  0  .  .  0 


o  B„ 


m=2 


9. 

0.  -0  .  .  B 
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m=M 


+i(2  (a/h)  (p/p  )  (c/c  ) 

_ § _ £_ 

/l-in 


m=l 

Bg  0  .  .  0 


0  B 
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0.  .0  .  .  B 
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nm=l 


nm=2 


nm=M 


nm^l 

?r  =2 
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nm=M 


(II.B-16) 


3.  The  Frame-Stiffened  Cylindrical  Shell 

The  effect  of  N^.  internal  stiffening  ring-shaped  frames  is  modeled 
here  as  a  set  of  initially  unknown,  radially  directed  (normal  to  the  shell 
surface)  modal  reaction  forces.  The  true  drive  may  be  taken  as  acting  on  the 
shell  through  a  particular  ring  frame,  or  as  acting  directly  on  the  shell  skin. 
The  procedure  for  the  solution  of  the  frame  reactions  requires  that  displace¬ 
ments  at  every  axial  frame  location  be  computed  due  to  each  virtual  load  in 
the  presence  of  fluid  loading.  An  opposite  system  of  virtual  reaction  forces 
acts  on  the  ring  frames,  for  which  a  second  set  of  equations  relating  frame 
displacements  to  reactions  may  be  written  with  knowledge  of  each  frame's 
structural  impedance.  Finally,  equating  shell  and  frame  displacements  at 
their  interface  yields  a  linear  system  of  equations  for  the  unknown  reactions. 

In  the  present  model  we  enforce  continuity  of  radial  and  tangential 
displacements  between  shell  and  frames  (w  and  q  in  Eq.  II.B-1);  so  that  for 
the  shell  not  only  must  we  calculate  the  radial  influence  of  the  superposed 
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system  of  virtual  radial  loads,  but  also  the  radial  influence  of  virtual 
tangential  forces  (equal  by  reciprocity  to  also  needed  tangential  dis¬ 
placements  due  to  radial  forces) ,  and  tangential  displacements  for  virtual 
tangential  forces.  However,  for  typical  submarine  shell  parameters  the 
latter  may  be  assumed  to  be  negligible  relative  to  the  others.  Therefore, 
if  is  used  to  denote  the  number  of  stiffening  frames,  the  order  of  the 
system  to  obtain  the  reaction  vector  is  2Nf;  and  each  of  the  2Nf  x  2Nf 
system  entries  is  a  combined  fluid-loaded  admittance  coefficient  for  the 
shell  and  frames.  The  known  right-hand  side  is  a  function  of  the  actual 
applied  drive.  After  the  radial  components  of  the  frame  reactions  have 
been  so  determined,  they  are  used  together  with  the  external  drive  (if 
applied  directly  to  shell  skin)  in  Eqs.  II.B-lc  as  the  total  system  of 
Fr(<f>/Z)  external  forces  driving  the  shell. 

We  calculate  the  modal  radial  velocity  and  surface  pressure  vectors 
for  the  cylindrical  part  of  the  boundary  by  solving  the  coupled  systems 
in  Eqs.  II.A-13,  14  and  II.B-15,  16  for  a  virtual  n-modal  load  acting 
in  frame  location  j.  Thus, 


TV 


nm=l 


nm=2 


,  V 


j 

nm-M 


1 


T 


( II . B-17a) 
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nm=l 


nm=2 


#  P 


j 

nm=M 


( II . B-17b) 


is  known  due  to 


fT>  \ 

<5(z-z.)  Z(  — - — J  cosnd) 

J  n \  pc  / 

where  the  r  subscript  on  F"1  denotes  the  radial  direction. 

n 

The  corresponding  tangential  modal  velocity  vector  Iv"1  ,  ,  v"1  „, 

nm=l  nm=2 

v3  ]  may  be  calculated  from  the  relationship  (v  =-iuq  ) 
nm=M  nm  nm 
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(a° '  \  ^ 


v  =1—  A 
nm  6 


v  +  .inUAI  {Fj  /a  _  pi  }  (II. 

nm  / , — r  r  nm 

p  c  v  1-in  n 
S  P 


so  that  radial  and  tangential  velocities  may  be  computed  at  axial  frame 
locations  z^ ,  i=l,  ...  ,  ,  through 


V  (z. ) ,  v  (z. )  =  £  (V  ,v  )  sin  (z.  +  L/2a) 

ni  ni  mnmnm  L  1 


(II. 


Thus,  defining  the  admittance  coefficients  c  c  by 

n. . ,  n , 


il  ID 


VJ (z. ) 
n  i 


(A 


(II. 


'th 


V  (z. ) 
n  l 


where  the  t  subscript  on  denotes 

the  total  velocity  at  a  she?!  surface 


tangential,  we  may  therefore  write 
point  z^  is  given  by 


V  (z.  ) 
n  l 


(II 


v  (z ,  ) 
n  l 


(II 
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In  the  above  it  is  assumed  that  the  applied  external  drive  acts 
indirectly  on  the  shell  through  a  frame  so  that  all  shell  forces  are 
virtual.  For  the  case  when  a  radial  drive  distribution  acts  directly  on 
the  shell  the  right  sides  of  Eqs.  II.B-21a,b  would  contain  respectively 
the  additional  terms 


rn  o 

-  c  . 

pc  ni 


o 

rn  -o 

-  c  . 

pc  ni 


where  c  .  denotes  the  loading  distribution's  radial  influence  at  position  i, 
ni 

and  c  .  the  tangential  influence, 
ni 

4.  Frame  Equations 


The  statement  of  Newton's  second  law  for  the  j  frame  takes 


the  form 


F:  /a 
r 

=  Z?.  V  (z.)+Z?_  V  (2.) 
pc  11  n  j  13  n  3 

n  n 


(II.B-22a) 


^t  /a 

— - =  zl  V  (z,)+z5j  v  (2.) 

pc  31n  n  3  33^  n  3 


(II . B-22b) 


where  the  frame  impedances  ,  ...  ,  have  been  derived  by  Harari  ,  and 

are  listed  here  in  Appendix  Il'?B(2).  The  net  force  F^_  is  given  by 


for  drive  on  frame 


for  drive  not  on  frame 


A  similar  expression  applies  for  F  .  Equations  II.B-22a,b  may  be  inverted 

in  order  to  obtain  explicit  forms  for  V  ,v  ,  which  then  may  be  set  equal  to 

n  n 

shell  displacements  to  yield  the  final  system  for  the  unknown  reactions  in 

^  ,  .  ...  comb 

terms  of  a  combined  admittance  matrix  [A]  : 
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As  earlier  stated,  iteration  schemes  such  as  those  proposed  here 

2 

require  also  a  number  of  operations  proportional  to  N  ,  the  precise  number 
2 

BN  depending  on  the  convergence  rate,  i.e.,  on  the  number  of  matrix-vector 

products  required  to  achieve  the  desired  level  of  accuracy  (each  matrix- 

2 

vector  product  is  an  N  operation).  On  the  other  hand,  exact  standard  solvers 

3 

such  as  Gaussian  elimination  require  CN  operations,  with  C=0(1)  . 

The  following  ratio  of  Total  Computational  Efforts (TCE)  therefore 
describes  the  savings  potential  of  iterative  methods  over  exact: 


(TCE) 

(TCE) 


I teration 


Exact 


2  2 
AN  +BN 

2  3 

AN  +CN 


(III. E-la) 


-  (-—)•-  asN*»  (III. E-lb) 

C  n 

where,  again, the  first  term  in  both  numerator  and  denominator  is  a  measure 
of  system  generation  effort,  and  the  second  terms  denote  actual  system-solving 
times.  Thus,  substantial  savings  could  result  wherever  the  needed  number  of 
iteration  steps  is  much  less  than  the  system  order N  . 
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where 


[A]  =  tcj]  -  [IS-bJ][1/C^]  (III.D-8) 

From  here  on  the  iteration  technique  follows  the  pattern  of  technique  2 

previously  discussed  for  n=0:  For  the  k=0  starting  step  we  set  x  (M+l)  = 

s  k=0 

x9(M+2)  =  0  and  obtain  {x^}  from  Eq.  III.D-7  by  solving  the  latter  in 

(inner)  Jacobi  steps  with  the  [B^] -matrix  corresponding  to  [A]  above.  The  cap 

mode  Eq.  III.D-5c  is  then  solved  either  exactly  or  approximately,  depending  for 

a  more  general  case  on  the  order  of  the  partitioned-of f  [Cl  subsystem.  For 

the  present  case,  [C]  is  2x2  and  an  exact  solution  is  used  in  Chapter  IV 

T  k=0 

applications.  With  these  values  of  [x  (M+l) ,  x  (M+l)]  on  the  RHS  of 

s  k=l  ^  ^ 

Eq.  III.D-7,  a  new  {x^}  vector  is  obtained  in  iteration  steps  and  so  on. 

As  in  Method  2  for  n=0,  the  whole  procedure  here  is  repeated  outer  steps. 

Also,  as  the  alternate  algorithm  for  n=0,  this  one  for  n>l  appears  to  be 
nonstationary  but  again  with  a  necessary  convergence  criterion:  S(B~)  <  1. 

This  second  algorithm  for  n>l  was  applied  successfully  whenever  Method  1  was 
found  to  diverge. 

E.  Economic  Advantage  of  the  Iterative  Methods 

The  total  cost  of  solving  fluid-structure  interaction  problems  has  two 
clearly  separate  contributions:  (1)  the  generation  of  the  systems,  i.e., 
filling  up  every  matrix  entry;  and  (2)  the  solution  process  applied  to  these 
systems.  In  a  finite-difference  or  finite-element  formulation  N  control 
points  are  chosen  to  represent  the  surface  of  the  structure  and  the  generated 
fluid  and  structural  systems  therefore  contain  NxN  elements.  The  work  needed 
to  fill  up  each  entry  depends  on  the  character  of  the  differential  equation 
describing  the  structure  and  on  the  numerical  integration  techniques  chosen 
to  perform  local  evaluations  of  the  Helmholtz  integral;  hopefully,  both  oper¬ 
ations  can  be  accomplished  with  many  fewer  than  N  multiplications  and  additions 

for  typically  large  systems.  Thus,  the  total  cost  of  system  generation  may  be 
2 

denoted  by  AN  ,  where  A  is  a  constant. 
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2 .  Method  2 

s 

Referring  back  to  Eq.  II.B-16,  we  recall  that  matrix  [C^]  ,  the 
structural  modal  admittance  matrix  for  modes  on  the  shell  surface  alone 
excluding  cap  modes,  is  diagonal  and  thus  readily  invertible  if  necessary. 

This  is  due  to  the  separability  of  this  part  of  the  geometry.  The  second 
iterative  algorithm  for  n>l  presented  next  exploits  this  feature  of  the  model 
formulation  without  invalidating  our  basic  policy  of  noninversion,  since  the 
latter  is  cost-effective  mostly,  and  thus  strictly  observed,  for  the  fluid 
system.  More  importantly,  since  structural  systems  resulting  from  the  in-vacuo 
finite-element  analysis  of  a  general  structure  normally  are  in  terms  of  an 
impedance  rather  than  an  admittance  matrix, no  inversion  is  necessary.  Also, 

g 

because  such  matrices  are  typically  heavily  banded,  the  product  [I-B^lfl/C^] 

that  arises  in  the  following  analysis  would  involve  many  fewer  operations  than 

3  s 

the  N  otherwise  required  if  the  impedance  matrix  [1/C^]  were  fully  populated. 

Following  an  approach  similar  to  that  for  n=0's  second  technique,  we 
write  that  now  for  n>l 

{ x® }  =  [C^Hxp  +  {bS}  (III  .D- 5a) 


[IS-B®){x®}  +  {(I-B  ).  }x  (M+l)  +  {  (I-B  )  .  mjJx,(M+2) 

2  2  2  l , M-f  1  2  2  i,M+2  2 

=  rc^HxJ} 


r  ~\ 

x2(M+l) 

i=l  (C2'M+l,i  xi(l)  “  (I_B2)M+l,i  *2(l) 

[C]  < 

►  =  < 

> 

M 

x  (M+2) 

.X  (C_)„  _  .  x  ( i )  -  ( I  —  B  )  _  .  xn(i) 

v  2  J 

1^1  =  1  2  M+2,l  1  2  M+2,i  2  J 

g 

Eq.  III.D-5a  now  is  inverted  to  solve  for  fx^): 


(III .D-5b) 


(III. D- 5c) 


{x®}  =  t i/c^] { xj}  -  r i/c®] { bs> 

which  upon  substitution  into  III.D-5b  yields 


(III.D-6) 


[A]  {xj}  =  -  [I^B^Hl/cjHb3}  +  {(I-B2).  m+1}x2(m+1) 


2  i ,M+2  2 


(III.D-7) 
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(III.C-13) 


Eq.  III.C-13  may  now  be  solved  iteratively  as  was  Eq.  III.C-11,  i.e., 
in  ^  iteration  steps.  The  procedure  is  stopped  when  k=K^  steps. 

*•'  g 

We  note  that  while  [B  ]  is  the  same  matrix  as  for  the  k=0  iteration, 

the  effective  dr^ve  vector  {b  }  is  now  clearly  different.  The  method 

therefore  is  non-stationary .  As  indicated  in  Ref.  29  ,  however,  the  fact  that 
-  s 

[B  ]  is  constant  vs.  k  keeps  the  convergence  analysis  for  this  subsystem  simple, 
and  a  necessary  condition  to  converge  now  is  that  S(BS)  <  1,  which  is  a 
different  criterion  from  that  obtained  previously  for  method  1  given  the  basic 
differences  in  system  matrix  partitioning  and  the  fact  that  the  total  scheme 
may  now  actually  be  semi-exact.  We  note  that  this  second  iteration  scheme  uses 
two  nested  algorithms:  The  outer  is  performed  times  and  corresponds  to  the 
k  counter;  the  inner  one  is  executed  K2  times  for  each  value  of  k. 

D.  Iteration  Algorithms  for  n  >  1 

Two  iteration  procedures  also  are  developed  for  the  systems  of 
Eqs.  III.B-2a  and  III.B-2b.  The  first  of  these,  however,  is  fundamentally  no 
different  from  the  method  1  for  n=0,  and  the  only  purpose  here  is  to  show 
briefly  how  the  analysis  for  convergence  still  follows  the  earlier  format 
though  the  pressure  and  velocity  vectors  now  differ  in  order,  and  also  one  of 

g 

the  coupling  matrices,  the  [<5  ]  matrix  of  Eq.  III.B-26,  is  nonsquare.  The 
second  algorithm  for  n>l  similarly  involves  some  of  the  same  elements  as 
method  2  for  n=0,  but  also  applies  a  novel  manipulation  not  previously  employed 
for  one  of  the  system  submatrices.  In  actual  applications,  this  second  algorithm 
was  observed  to  converge  when  the  first  one  did  not. 

Following  earlier  nomenclature,  we  now  write  Eqs.  III.B-2a  and  III. B- 2b  in 
equivalent  form  suitable  for  iteration: 

{xl}  =  [Ci1{x2}  {h)S}  (III.D-la) 

fx2}  =  [B2]{x2}  +  (C2lfx®} 
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X  T 

To  start,  fxj_'x2^  i-s  set  to  zero  so  that  the  k=0  iterate  of  [q^,q^]  , 

ooT  ooT  sT 

denoted  [g^g^  ,  becomestg^^^]  =  [b  ,01  in  Eq.  III.C-8a,  which  then  reads 
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■Tv  — kJ 

r  s"> 

i  -c. 

g„ 

b 

i 

1 

< 

►  =  < 

►  =  < 

l 

s  s  s 

s 

-C .  I  -B 

0 

2  2 

v  2  J 

k_  J 

(III.C-10) 


Eq.  III.C-10  may  then  be  solved  by  Jacobi's  method  (or  by  other,  more 
sophisticated  techniques)  in  a  set  of  J-iteration  steps.  We  write, 


k=0 


r  s 

xi 


k=0 


=  [BS]  < 


+  {bS} 


k=0 


(III.C-11) 


^2J 


j-1 


where 


r 


-s,  -1 


[B1  =  [ D  ]  <  ( D  ]  - 


-c: 


**-L 


(III .C-12a) 


9  " 


{bS}  =  [DS]  1 


(III.C-12b) 


and  where  [D  ]  is  a  diagonal  matrix  containing  Eq.  III.C-10' s  diagonal  elements. 

Again,  the  iteration  in  Eq.  III.C-11  is  carried  out  steps,  after  which  the 

s  s  Tk=0  m  ^  Tk=0 

results  so  obtained  for  fx1»x2l  are  used  to  computed  [glfg2l  from 

Eqs.  III.C-9e  and  III.C-9f  for  use  in  Eqs.  III.C-8b.  The  latter  may  now  be 

s  s  Tk=0 

solved  iteratively  in  the  same  fashion  as  was  [k^jX^]  ,  or,  if  this  subsystem 

is  conveniently  chosen  to  be  small  enough  as  in  the  present  case,  where  it 

~  ~  Tk=0 

contains  only  the  two  cap  unknowns,  an  exact  solution  for  [x  ,x2)  may  be 

«  -  ~  Tk=0 

just  as  economical.  These  values  of  [x.,x  ]  are  then  used  to  calculate, 

using  Eqs.  III.C-9c  and  III.C-9d,  new  values  of  g^,g2,  denoted  g^  ,  g2  for 

another  k-iteration  with  Eqs.  III.C-8a,  which  then  reads: 


j.;-.  ^ 
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where 


x1(M+l) 


■'v 


{x^  = 


x1(M+2) 


x2(M+1) 


{x2}  = 


x  2  ( M+  2 ) 


<*!>  =  - 


{ai,M+l}xl(M+1)  "  ^ai,M+2}xl (M+2)  +  { (C15 i,M+lix2(M+1) 
+  {(Cl>i,M+2}x2(M+2)  +  {bS} 


tg2)  = 


{(I-B2)i,M+l}x2(M+1)  -  {(I"B2)i,M+2}x2(M+2) 


+  {(C2>i,M+l}xl(M+1)  +  {(C2)ifM*2,xl(M+2) 


i=l  ^C1 ^  M+l , i  X2(l) 


{g, }  = 


M 


^i=l  (Cl5M+2,i  X2(l). 


M 


i=l  {I~B2)M+l,i  X2(l) 


M 


{i2}=- 


I  M 

U 


ill 


.  x  (i) 

2  M+2 , l  2  j 


i  =  l  ^2)M+2,iXl(l) 


(III .C-8a) 


(III.C-8b) 


(III .C-9a) 


(III .C-9b) 


(III.C-9c) 


(III .C-9d) 


(III.C-9e) 


(III .C-9f ) 
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M 


[a]  J 


"x  (M+l)"' 
1 

"x2(M+lT 

^II  (Cl)M+l,i  X2(i)' 

< 

>  -  ICX1  < 

►  +  < 

> 

M 

x  (M+2) 

x  (M+2) 

. £_  (C. )  .  x^(i) 

L  1  J 

l  2  J 

^i=l  1  M+2,i  2  J 

[IS-B^]{xS2}  +  {d-B2)ifM+1}  x2  (M+l )  +  {d-B2)ifM+2}  x2  (M+2) 


=  icjH*;}  +  {(c2)i>M+1>  xi(m+d  +  {(c2)i(M+2}  *x 


(M+2) 


[T-b2] 


f x  (M+l ) > 

r.T  (I-Bl  ,  .  X®(i) 

2 

1=1  2  M+l , i  2 

►  =  -  < 

) 

M  ,  [ 

x  (M+2) 

.E,  (I-B.)„  „  .  x®(i) 

L 2  J 

^1=1  2  M+2,i  2  J 

rx1(M+l)  ^ 


M 


+  [c2]  < 


i=l  (C2)M+l,i  xl(l) 


,Xl(M+2) J  U  <C2) M+2 , i  Xl(i)  J 


(III. C- 7b) 


(III.C-7c) 


(III .C-7d) 


where  {a.  „  {a.  .}  are  M-element  vectors  made  up  of  the  first  M  entries  of 

1 , M+l  i ,M+2 

the  [a]  matrix's  (M+l)tl:l  column  (Eq.  III.B-la).  Similar  definitions  apply  to 

{ (C  ) .  },  { (I-B„) .  },  etc.  Eqs.  III.C-7a  and  III. C-7c  have  order  M. 

1  i,M+l  2  i,M+2 

s  s 

Also,  it  will  be  recalled  that  the  upper  left  MxM  submatrix  [I  -B  ]  of  [a]  is 

the  identity  matrix  since  [B®]  =  0.  Systems  III.C-7b  and  III.C-7d  have,  for  the 

1  M  s 

present  problem,  an  order  of  2.  A  summation  such  as  (I-B_)„  ,  .  x„(i) 

1=1  2  M+l ,i  2 

t  ri  g 

denotes  the  inner  product  of  the  (M+l)  row  of  matrix  [I-B2]  with  vector  (x2). 

The  2x2  matrices  ci  ,  J>b2 ,6 , C2  contain  the  bottom  rightmost  four  elements(M+l ,M+1 ) , 
(M+l, M2),  (M+2, M+l),  (M+2, M+2)  of  matrices  a it  is  important  to 

point  out  that  such  partitioning  of  the  two  original  systems  should  be 
mathematically  completely  arbitrary.  For  the  physical  problem  at  hand,  however, 
such  splitting  sets  aside  the  M-subvector  of  shell-surface  effects  from  the 
two-element  cap  subvector  by  placing  these  quantities  on  opposite  sides  of 
the  equations.  If  more  cap  modes  had  been  used  in  the  model,  the  two  smaller 
subsystems  would  have  been  augmented  accordingly.  Eqs.  III.C-7a  through  III.C-7c 
may  also  be  written  in  combined  form: 
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one  realizes  that  expansions  III.C-2  through  III.C-4  are  also  obtained,  to  the 
identical  order  of  error  indicated,  by  simple  Jacobi  iteration  of  the  combined 
system: 


(III.C-5) 


~~k-l  c 
x  =  Bx  +  b 


(III.C-6) 


where  x  is  now  a  2 (M+l) -element  vector,  etc.  Thus  identifying  the  equivalence 
of  algorithm  III.C-1  to  that  corresponding  to  Eq.  III.C-5  establishes  its 
convergence  criterion,  which  is  now  clearly  that  S(B)<1.  In  summary,  fixed- 
point  iteration  of  joint  structure- fluid  systems  converges  if  the  matrix  B, 
composed  of  the  "self"  matrices  B^,B2  along  its  diagonal  blocks,  and  of  the 
"coupling"  matrices  C^,C  as  off-diagonals,  has  a  spectrum  all  less  than  one 
in  magnitude.  Consistency  proof  of  algorithm  III.C-1  (i.e.,  that  if  for  some 


.  *  k*  k* 

k=k  ,  xx  ,x2 


x*  =  (x  ,x  )exact,  then  =  xfxact  also)  follows  trivially  from 

Z  X  Z  X  $  Z  1/4 


the  standard  single-system  proof.  Similarly,  although  at  first  algorithm 

III.C-1  does  not  appear  to  be  stationary,  i.e.,  that  for  both  x1  and  x2  the 

effective  "self"  matrix  and  associated  known  RHS  appear  to  be 

functions  of  the  iteration  counter  k  rather  than  being  k- independent  quantities, 

the  equivalent  system  of  Eqs.  III.C-5  and  III.C-6  in  fact  shows  that  the 

algorithm  is  stationary  and  that  for  the  combined  system  B  „wb  „  are  just 

™  eff  eff 

— 

the  constants  B,  [b,o]  ,  where  T  denotes  transpose.  As  previously  stated,  once 
the  iteration  scheme  has  been  recast  in  terms  of  Eq.  III.C-5,  the  latter  may  be 
changed  to  other  techniques  with  accelerated  convergence,  e.g.,  SOR. 

2 .  Method  2 

Eqs.  III.B-3a  and  Til. B- 3b  may  also  be  written  as 


-n! 


CIS-B®|{x®>  +  M+1}xl(M+l)  +  («ifM+2}x1(M+2) 

=  [C^HxJ}  +  {(C1)ifM+1}x2(M+1)  +  {(Cl)i>M+2}x2(M+2)  +  {bS} 


(III .C-7a) 
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C.  Iteration  Algorithms  for  n=0 
1.  Method  1 

To  simplify  the  nomenclature,  we  drop,  in  some  of  what  follows 
next,  the  curly  and  square  brackets  from  quantities  defined  in  Eqs.  III.B-3, 

III.B-4  and  III.B-5.  The  first  algorithm  we  wish  to  explore  for  Eqs.  III.B-3a  and 
III .B-3b,  is 

k  k-1  k-1 

x^  =  B^x^  +  +  b  (Ill.C-la) 

k  k-1  k-1 

x2  =  B2x2  +  C2Xl  (Ill.C-lb) 


To  begin,  one  sets  x°  =  0  and  =  b.  The  solution  x°  =  b  does  not  really 

-1 

correspond  to  the  in-vacuo  response  of  the  structure  because  b  is  not  a  f,  but 


.-1, 


instead  only  f.  The  first  few  iterations  go  as  follows: 


k  =  1 


k  =  2 


k  =  3 


X1  =  (I+Bi)b 


x"  =  C  b 
2  2 


*1  -  d^BiW^b 


*2  -  (C24C2Bl+B2C2,b 


*1  '  <WVBXlb  +  <ClC2+BlClC2+ClB2C2+Cl°2B1,b 


*2  '  <C2tC2Bl+B2C2+B2C2+B2C2Bl+C2BJ+C2ClC2)b 


(III .C-2a) 

(III. C- 2b) 

(III .C-3a) 

(III .C-3b) 

(III .C~4a) 

(III .C-4b) 


At  first  the  analysis  for  convergence  of  x^  and  x2  given  the  expansions 
Eqs.  III.C-2  through  III.C-4  appears  complicated  due  to  their  coupling.  For 
example,  although  the  B^b  part  of  the  solution  for  x^  is  reminiscent  of 
iteration  by  Jacobi's  method,  there  is  also  the  C^C^b  +  •••  Part  which 
necessarily  changes  the  criterion  of  convergence  from  S(B  )<1  to  something  else 
(we  follow  the  standard  nomenclature  of  Ref.  29,  where  8(8^)  stands  for  the 
spectrum  of  matrix  B^;  thus  S(B^)<1  states  that  all  eigenvalues  of  B^  must  be 
less  than  1  in  magnitude) .  This  situation  is  greatly  simplified,  however,  when 
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purposes  of  iteration,  Eqs.  III.B-la  and  III.B-Ib  for  n=0  are  now  written  as 
follows: 


{x.^} 

=  [B1]{x1>  +  [C1]{x2)  +  {b} 

(III ,B-3a) 

{x2} 

=  [b2Hx2}  +  [C2]{Xl} 

(III. B- 3b) 

[B1] 

=  [D  ] -1 [D  -a] 

a  a 

(III .B-4a) 

[C1] 

=  [d  r1!?,] 

a  1 

(III .B-4b) 

{b}  = 

=  to  ]_1[e„Hf} 
a  2 

(III.B-4c) 

[B2] 

=  [D  ] _1 [D  -y] 

Y  y 

(III.B-4d) 

[C2] 

=  [Dy]-1[6] 

(III.B-4e) 

where  stands  for  the  diagonal  matrix  composed  of  the  diagonal  elements  of 
[a] ,  and  [D^]  for  that  corresponding  to  [y] . 

For  n>l,  one  may  similarly  write 


{ x, }  =  [C, ] {x„}  +  {b  } 


(III. B- 5a) 


{x2}  =  [B2Hx2)  +  [62]{xJ}  (III  ,B-5b) 

where  {bS}  =  {fS},  and  [B  ]  is  still  given  by  Eq.  III.B~4d.  The  (M+2)xM  matrix 

^  —Is 

[C2]  is  given  by  the  product  [D^]  (<5  ]. 

For  both  the  n=0  and  n>l  formulations  we  now  develop  two  separate  iteration 
algorithms  and  discuss  their  convergence.  For  example,  for  the  n=0  problem  we 
show  that  the  two  algorithms  place  upon  the  joint  system  different  analytically 
derived  convergence  criteria,  so  that  if  one  algorithm  fails  the  other  may 
possibly  converge.  The  n>l  problem  is  treated  similarly.  We  will  show  in 
Chapter  IV  that  over  the  frequency  range  explored  at  least  one  of  the  iteration 
techniques  always  converged  to  the  exact  solution. 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


which  contains  a  free  convergence-optimization  parameter  not  analytically 
available  in  the  J-method.  However,  no  attempt  is  made  here  to  adapt  such  SOR- 
type  techniques  to  the  basic  algorithms,  and  we  merely  note  that  the  savings 
realized  as  reported  in  Chapter  IV,  while  significant,  could  thereby  be  made 
even  greater. 

B.  System  Equations 

Equations  II.A-13  and  II.A-14  relate  surface  modes  of  fluid  pressure 
to  surface  modes  of  fluid  particle  velocity  for  the  nonseparable  finite-cylinder 
geometry.  Eqs.  II.B-15  and  II.B-16  relate, for  the  structure, vibrational 
response  to  an  applied  drive  in  the  presence  of  the  surrounding  fluid,  which 
contributes  an  unknown  load  distributed  over  the  axial  surface  of  the  shell 
and  end  caps.  Equating  the  surface  response  vector  to  the  surface  fluid 
particle  velocity  yields  the  final  coupled  systems  from  which  {v}  and  {p}  may 
be  solved  for  a  given  drive  load  configuration  {f}.  Eqs.  II.A-13  and  II.B-15 
apply  for  circumferential  n-mode  n=0,  and  Eqs.  II.A-14  and  II.B-16  apply  for 
n>l.  Thus  in  summary,  we  have  that  for  n=0. 


[a] {v}  =  te1i {P}  +  ce2i {f> 

[y]{p>  =  [<5Hv} 


(III.B-la) 


(III.B-lb) 


where  every  matrix  contains  (M+2)x(M+2)  elements  (each  row  and  column  contains 
M  shell  surface  models  plus  one  cap  mode  for  each  end  of  the  cylinder,  and  each 
vector  contains  M+2  elements.  For  n>l,  we  have 


{vs}  =  res]{Ps}  +  [8®}{f s) 


(III.B-2a) 


tyHp}  =  C 6 S 1  { v  } 


(III .B-2b) 


where  [y]  and  {p}  still  have  (M+2)x(M+2),  and  M+2  elements,  respectively,  but 
where  now  {vS}  denotes  the  shortened  M-element  modal  velocity  vector  which 
contains  only  shell-surface  modes;  we  recall  tha+  this  is  the  result  of 
assuming  that  for  n>l  the  caps  are  rigid.  Similarly,  the  s-superscript  on 
B2,  and  f  indicate  a  rank  of  M  rather  than  M+2.  Matrix  [6S]  contains  M+2  rows, 
one  for  each  of  the  M+2  pressure  modes,  for  which  a  "diffracted"  cap  mode 
exists;  and  it  contains  M  columns,  one  for  the  effect  of  each  of  the  M  shell- 
surface  velocity  modes  on  each  pressure  mode,  including  cap  pressures.  For 


L 
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III.  ITERATION  TECHNIQUES  FOR  COUPLED  STRUCTURAL-FLUID  SYSTEMS  FOR  A  FINITE 
CYLINDER 

A.  Introduction 

The  main  purpose  of  the  present  study  is  to  show  whether  fully 

populated  structural-fluid  systems  as  would  arise  from  a  finite-element  or 

modal  analysis  of  a  nonseparable  geometry  may  be  solved  iteratively  and  thus 

more  economically  than  by  the  standard  exact  method,  which  effectively  requires 

the  inversion  of  at  least  one  large  matrix.  For  N-order  coupled  systems  the 

3 

exact  solution  requires  a  number  of  operations  proportional  to  N  ,  while  the 

2 

corresponding  number  for  iterative  techniques  is  only  N  .  In  our  earlier 
investigation  for  the  spherical  shell1  we  demonstrated  that  well-established 
iterative  techniques,  e.g.,  Richardson's  method  and  the  so-called  method  of 
summability  (Ref.  29,  p.  345),  could  be  applied  to  solve  those  versions  of  the 
Helmholtz  and  dynamic  response  equations.  However,  due  to  the  separability  of 
the  spherical  geometry,  all  matrices  were  effectively  diagonal  and  the  viability 
of  this  promising  approach  went  untested  for  the  more  fundamentally  complex 
general  case. 

In  this  chapter  we  develop  iterative  schemes  to  solve  the  finite  cylinder 
problem  as  a  representative  of  a  wider  class  containing  off-diagonal  elements. 

For  each  procedure  we  discuss  qualitatively  conditions  for  successful  iteration 
and  attempt  to  explain  why  a  specific  technique  may  work  when  another  fails. 
Actual  numerical  results  and  comparison  to  exact  solutions  are  deferred  to 
Chapter  IV,  where  we  also  tabulate  computational  times  vs.  number  and  type  of 
iterations.  We  find  that  the  stability  of  computations  often  hinges  on 
judicious  system  subpartitioning,  a  concept  which  has  been  commonly  applied, 

though  in  somewhat  different  context,  to  the  otherwise  difficult  treatment  of 

t  .  ...  ,  .  37,38,39 

certain  ill-conditioned  matrices. 

One  general  comment  applies  to  the  convergence  rate  of  essentially  all 
the  algorithms  developed.  Initially  we  show  that  the  resulting  expansion  for 
parts  of  the  iteration  algorithms  may  be  interpreted,  after  proper  partitioning 
of  matrices  and  vectors,  as  an  application  of  fixed-point,  or  Jacobi,  iteration 
(Ref.  29,  p.  71).  In  some  cases,  the  convergence  rate  of  such  J-method  solutions 
could  be  improved  substantially  by  a  recasting  to  yield  equivalent  expansions 
corresponding  to  the  Successive-Over-Relaxation  (SOR)  method  (Ref.  29,  p.  73), 
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N  N 

„  comb  r  1  2  f  1  2  „  f,T 

[A]  [F  ,  F  ;  Ffc  ,  F  ,  . . . ,  Ffc  ] 

n  n  n  n  n  n 

^  ^  r 

/  [0,0, . . ,Ai*,0. . . ,0;0, . . . . . ,0] 


[-c, » -c  , . . , -c  ;  -c, ,  -c  ,  ...,-c  J 


1  2 


1  2 


(II.B-23) 


where 


z1  Z  -Z1  Z1 
11  33  13  31 

n  n  n  n 


B.  =  — ; - ; — ; - — 

1  7 ^  Z1  — Z1 

11  Z33  Z13  31 

n  n  n  n 


(II.B-24) 


(II.B-25) 


The  right  sides  of  Eq.  II.B-23  are,  (1)  radial  drive  applied  to  the  i  frame, 
(2)  radial  drive  acting  directly  on  the  shell.  In  the  former  case  the  model 
assumes  that  the  applied  drive  distribution  is  concentrated  at  an  axial  lo¬ 
cation,  while  in  the  latter  the  radial  drive  has  an  arbitrary  axial  distri¬ 
bution.  For  both  cases  the  circumferential  drive  distribution  is  arbitrary. 

Again,  upon  solution  of  Eq.  II.B-23,  the  combined  radial  loading  system 
on  the  shell  (reactions  plus  drive)  becomes  the  effective  shell  drive  which 
together  with  fluid  loading  appears  forcing  the  structural  systems  of 
Eqs.  II.B-15  and  16). 
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IV.  APPLICATION  OF  THEORY 

A.  Introduction 

The  modal  solution  technique  developed  in  Chapter  II  is  now 
applied  to  the  idealization  of  an  actual  structure  for  which  extensive 
far-field  measurements  exist.  Comparisons  are  made  between  theory  and 
experiment.  The  iterative  techniques  of  Chapter  III  are  also  applied 
and  their  results  are  compared  to  the  exact  solution. 

The  cylindrical  structure  shown  in  Figures  IV.la,b  is  one  of  five 
equal  segments  spliced  to  make  up  the  total  physical  model  in  Ref.  40, 
which  is  a  recent  compilation  of  related  experimental  acoustic  and  re¬ 
sponse  data  collected  by  the  Penn  State  group  of  Burroughs,  Hayek, 

Hallander,  and  Bostian,  in  coordination  with  the  ONR.  Figure  IV.  2 
shows  the  simplified  structure  used  for  our  calculations.  Only  the 
large  joint-band  frames  have  been  kept  as  structural  discontinuities, 
while  all  non- axi symmetric  details,  including  axial  stringers,  have  been 
eliminated.  The  simplified  structure  models  what  Ref.  40  calls  the  single, 
or  inner,  shell  configuration.  The  predictions  reported  here  are  for  a 
radial  point  drive  applied  either  directly  to  the  shell  skin  at  the  cylinder's 
midsection  point,  or  at  on:  jf  the  two  frames  nearest  the  center.  Structural 
parameters  corresponding  to  Figure  IV. 2  are  listed  in  Table  IV. 1. 

For  reasons  outlined  below  having  to  do  with  system  size  and  a  limited 
computational  budget,  calculations  were  not  carried  out  by  the  present  method 
beyond  500  Hz.  In  order  to  continue  the  comparison  of  measurement  to  pre¬ 
diction  to  2500  Hz,  the  fluid  model  was  replaced  by  that  for  an  infinite 
cylinder. 

For  clarity  of  presentation,  the  discussion  and  physical  interpretation 
of  predictions  are  given  first  at  some  length  using  the  results  obtained  from 
solving  the  coupled  systems  exactly.  We  demonstrate  subsequently  that  in 
fact  all  such  exact  results  are  also  obtainable  by  iteration;  thus,  they  may 
be  viewed  as  converged  solutions.  The  calculations  describing  just  the 
iteration  process  are  then  displayed,  both  in  frequency-sweep  form  and  for 
directivity  patterns.  In  order  to  present  as  many  results  as  possible  with 
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a  minimum  of  repetition,  some  of  the  directivity  patterns  for  the  iterated 
solution  correspond  to  frequencies  which  are  different  from  those  shown 
under  the  section  for  the  exact  results.  We  conclude  that  the  iterative 
solution  of  fluid-loaded  structural  systems  is  both  feasible  and  cost 
effective,  at  least  in  the  low-frequency  range  where  fluid  loading  has 
a  global  effect  and  is  therefore  less  amenable  to  local  modeling. 

Since,  as  pointed  out  in  Chapter  III,  convergence  rates  do  not  de¬ 
pend  on  the  coupled  systems'  forcing  function,  iterated  solutions  will  be 
presented  only  for  the  drive-on-skin  case.  Iterated  solutions  for  the 
driven  frame  may  be  assumed  to  converge  similarly. 

B.  System  Size  vs.  Frequency 

The  modal  theory  of  Chapter  II  requires  that  the  series  of 
simply- supported  axial  modes  be  truncated  to  yield  coupled  systems  of 
finite  rank.  Here  we  determine  quantitatively  the  minimum  needed  number  M 
of  such  modes  in  terms  of  the  frequency-dependent  flexural  wavenumber  kf 
for  the  shell  plating.  The  minimum  system  size  M  must  capture  the  shell 
flexural  wavenumber  kf  and  if  possible  overshoot  it  by  a  certain  margin. 

The  inclusion  of  kf  in  the  m  sum  allows  a  fair  representation  of  the 
structural  far- field  effect  of  each  frame  on  the  others,  as  required  by 
off-diagonal  frame  admittance  coefficients  in  Eq.  II.B-20;  while  the 
additional  modes  are  needed  to  calculate  the  diagonal,  self  admittance 
elements.  Thus,  we  have  that 

M  >  —  •  k  a  (IV.B-1) 

7T  f 

With  Eq.  T.V.B-1  serving  as  the  criterion  for  system  size,  the  minimum 
matrix  order  M  vs  frequency  is  given  on  the  second  column  of  Table  IV. 2. 

For  reasons  of  economy,  only  the  range  20  Hz<f<500  Hz  was  investigated 
here  and  the  value  M  was  kept  fixed  at  30  for  all  calculations. 

Values  of  the  normalized  acoustic  wavenumber  ka  are  given  on  the 
third  column  of  Table  IV. 2.  We  note  that  for  every  frequency  the  below- 
coincidence  condition  k^a<ka  is  satisfied,  so  that  every  plating  wavenumber 
listed  has  a  subsonic  trace  speed  and  thus  does  not  contribute  directly  to 
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the  far-field.  The  coincidence  frequency  for  the  shell  skin  is  approx¬ 
imately  35.3  kHz  and  therefore  well  above  the  range  of  interest  here. 
However,  for  the  shell  in  question  the  axial  modal  wavenumber  mir/(L/a) 
is  given  roughly  by  m-(.14),  so  that  for  example  500  Hz  contains  m=l,2,3 
modes  which  are  supersonic  and  therefore  radiative;  the  rest,  4<m<30=M, 
are  subsonic.  The  relative  contribution  of  subsonic  and  supersonic  modes 
to  the  acoustic  field  has  been  discussed  qualitatively  in  some  detail  in 
Chapter  II,  Section  A. 

C.  Choice  of  Number  of  Circumferential  Modes 


In  the  theoretical  model,  stiffening  frames  represent  axial 
structural  discontinuities;  however,  since  no  circumferential  discontin¬ 
uities  which  would  render  the  structure  nonaxi symmetric  are  considered 
here  either  for  the  framing  or  for  the  skin,  a  structural  criterion  ana¬ 
logous  to  that  of  Eq.  IV. 1  need  not  be  imposed  on  the  circumferential 
mode  n.  The  only  requirement  on  the  number  of  n  modes  is  that  they  in¬ 
clude  all  those  which  radiate  efficiently,  i.e.,  all  n  satisfying  n<ka. 
Here  we  have  used  the  n=0,l  and  2  modes,  of  which  according  to  Table  IV. 2 
column  3,  only  the  first  is  supersonic  over  20  Hz  <f<500  Hz.  However, 
strong  n=0  response  is  not  anticipated  for  the  low  frequencies  considered 
and  all  three  modes  may  in  some  cases  contribute  comparably  to  the  far- 
field  signal.  Acoustical  radiation  from  subsonic  and  supersonic  circum¬ 
ferential  modes  was  also  discussed  qualitatively  in  Chapter  II,  Section  A. 


The  Effect  of  Stiffening  Frames  on  Acoustic  Radiation 


In  the  low  frequency  limit  frame  structural  impedances  become 
inertial  and  thus  vanishingly  small.  Internal  frame-skin  reactions  are 
similarly  zero  as  the  entire  structure  moves  in  its  fundamental  beam  mode 
As  the  frequency  increases,  frame-shell  impedance  mismatches  grow  and  so 
must  reaction  forces.  Frame  reactions  may  be  expected  to  be  greater  for 
those  frames  positioned  away  from  the  nodes  of  an  excited  m  shell  mode; 
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e.g.,  an  m=l  resonant  mode  will  generally  cause  greater  reactions  for 
the  two  frames  near  the  shell  middle  than  for  those  near  the  two  ends. 
Conversely,  the  m=5,  10,15. . .resonances  should  have  associated  smaller 
reactions  since  the  four  frame  positions  coincide  naturally  with  modal  nodes. 
Given  sufficient  structural  damping,  at  higher  frequencies  frame  reactions 
become  acoustic  points  (or  rings)  of  origin  wherever  the  shell-plating 
structural  far  field  they  excite  is  below  coincidence. 

E.  Exact  Solutions:  Numerical  Results 

E. 1  Summary 

All  far-field  predictions  presented  here  use  the 

spherical  coordinate  system  of  Figure  II. 1.  The  field  position  defined 

by  <j>=0,  iJj=90°  (beam  aspect)  is  in  line  with  the  drive  direction;  <J>=180° , 
o 

tj;=90  points  in  the  direction  opposite  the  drive.  Preliminary  spot 
checks  on  the  shape  of  the  acoustic  directivity  revealed  that  for  most 

frequencies  in  the  range  investigated  pressure  maxima  occurred  within 

o  o  o 

elevation  angles  i|i=65  and  85  ,  typically  in  the  vicinity  of  t|>=80  . 

Far-field  results  are  presented  at  20°  intervals  for  both  circumferential 

and  elevation  angles.  Since  directivity  patterns  were  fairly  smooth,  the 

levels  shown  for  i|)=80°  are  in  essence  equal  to  those  for  i^=90°,  which  are 

not  shown.  The  dimensional  distance  Ra  from  the  cylinder's  center  to  the 

far-field  point  is  45  ft;  thus,  k'Ra  varies  in  value  from  about  1  at  20  Hz 

to  just  over  25  at  500  Hz.  The  pressure  has  been  normalized  by  the  quantity 
2 

F/a  ,  where  F  is  the  magnitude  of  the  applied  drive. 

E. 2  Exact  Solution:  Frequency  Sweeps 

Figure  IV.  3  shows  frequency  sweeps  at  i(i=80O,  <(>=0  and  at 

o  o 

i|j=80  ,  c(>=180  for  the  drive  on  the  shell.  Figure  IV. 4  shows  similar  results 
for  the  drive  on  the  frame.  The  acoustic  field  for  the  drive  acting  directly 
on  the  fluid  is  given  in  both  by  the  force-dipole  curve.  The  latter  far 
field  is  antisymmetric  with  respect  to  the  (j>=90°  plane  and  its  levels  for 
d>=0  and  180°  are  identical.  In  Figure  IV.  3  the  spectrum  for  <{>=0  generally 
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stays  above  the  dipole  curve,  while  that  for  (j>=180°  seems  to  oscillate 

o 

about  it.  For  lower  frequencies  the  <j>=0,  180  curves  merge  in  agreemept 
with  the  expectation  that  the  structural  response  should  then  be  primarily 
rigid  body  with  an  associated  dipole-like  far  field.  In  the  extreme  low- 
frequency  range  bounded  above  by  the  fluid-loaded  fundamental,  the  system's 
total  impedance  is  stiffness-controlled  and  some  departure  from  the  dipole 
curve  should  be  expected.  This,  however,  is  an  artifact  of  the  model  since 
the  subject  actually  was  a  freely  suspended  hull. 

Figure  IV. 4  corresponds  to  the  driven-frame  case.  The  relatively 
small  differences  between  the  <P=0,  180°  curves  indicate  rigid-body  motion 
over  most  of  the  spectrum. 


Exact  Solution:  Directivity  Patterns  for  Drive  on  Shell 


Figures  IV. 5-8,  corresponding  to  the  drive-on-shell  case, 
show  predicted  directivity  patterns  for  a  few  of  the  peaks  in  Figure  IV. 3, 
<J>=0  spectrum.  Each  figure  shows  three  curves  corresponding  to  three 
different  elevation  angles;  each  curve  shows  far-field  pressure  vs.  circum¬ 
ferential  angle  ^  on  a  plane  normal  to  the  cylinder  axis.  Thus  for  example, 
the  <fi=0  value  of  the  \p=80°  curve  on  Figure  IV.  5,  roughly  -64  dB,  corresponds 

to  the  local  peak  at  100  Hz  in  Figure  IV.  3.  Similarly,  the  <J>=180° value  for 
o  o 

ijj=80  corresponds  to  that  for  <J>=180  at  100  Hz  in  Figure  IV.  3.  The  low 
values  obtained  near  c()=90O  (direction  perpendicular  to  the  drive  at  every 
elevation  angle) further  confirms  the  dipole-like  behavior  of  the  acoustic 
field.  Since  the  drive  acts  at  the  cylinder's  exact  center,  the  directivity 

patterns  of  Figures  IV.  5-8  are  symmetric  about  the  t|*=90°  plane;  e.g.,  the 

o  o  o 

pattern  for  t|)=100  is  the  same  as  that  for  ip=80  ,  and  that  for  ij>=140  is 

equal  to  that  for  iji=40°. 

The  gap  at  200  Hz  between  the  <p=0,  180°  curves  has  been  attributed 
to  n=0,2  contributions  comparable  in  magnitude  to  that  of  n=l.  Figure  IV. 6 
shows  the  200  Hz  far-field  pattern.  We  note  that  while  the  largest  pressure 
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still  belongs  to  the  ij/=80°  elevation,  other  elevations  may  dominate  at 
other  <p  positions. 

Figure  IV. 7  and  8,  respectively  for  300  and  400  Hz,  continue  the 
trend  towards  increased  complexity  with  increasing  frequency,  especially 
for  lower  elevations. 

Figure  IV. 9  shows  the  directivity  pattern  on  the  plane  containing 
the  cylinder  axis  and  the  (j>=0  point.  The  circumferential  angle  is  fixed 
at  (J>=0  and  the  elevation  angle  varies  from  the  position  tp=0,  which  looking 
to  the  cylinder  sees  only  the  right  cap,  to  ^=180°  which  sees  only  the 
left  cap.  Curves  are  shown  for  the  four  peak  frequencies  displayed  in 
Figures  IV. 5-8.  As  previously  stated,  all  curves  are  symmetric  with  re¬ 
spect  to  \p=90° .  For  each  curve,  the  ip=Q,  180°  values  denote  predicted 
radiation  levels  at  the  two  far-field  points  which  coincide  with  the 
cylinder  axis.  These  levels  increase  with  frequency  (as  expected) ,  though 
apparently  not  monotonically . 

Since  the  theoretical  model  allows  cap  motion  only  for  n=0,  low 
levels  at  ij>=0,  180°  for  100  Hz  are  consistent  with  the  earlier  assertion 
that  the  structure  primarily  responds  then  as  a  rigid  beam;  for  n=l  the 
caps  are  assumed  rigid,  and  low  but  nonzero  levels  at  iji=0,  180°  are  due 
then  to  diffraction  of  signals  emanating  from  the  shell  surface.  The 
higher  levels  observed  for  the  other  curves  point  indirectly  to  the  greater 
participation  of  the  breathing  mode  as  frequency  increases.  Finally,  as 
previously  mentioned,  the  highest  absolute  levels  were  predicted  near 
if/=80°  for  most  frequencies. 

E . 4  Exact  Solution;  Directivity  Patterns  for  Drive  on  Frame 

o 

The  peaks  in  the  calculated  spectrum  at  i^=80  shown  in 
Figure  IV. 4  for  the  drive-on- frame  case  correspond  approximately  to  fre¬ 
quencies  100,  200,  350,  and  500  Hz.  For  these,  directivity  patterns  are 
shown  in  Figures  TV.  10-14.  Figures  IV.  10-13  show  patterns  on  planes 
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normal  to  the  cylinder  axis;  Figure  IV.  14  shows  patterns  on  the  plane 

defined  by  the  cylinder  axis  and  the  point  cf>=0.  Unlike  the  drive-on- 

shell  (middle)  case,  these  radiation  patterns  are  not  symmetric  relative 

to  the  ij)=90°  plane.  Thus,  to  provide  a  better  overall  picture  of  the 

acoustic  field,  the  elevations  angles  displayed  have  been  changed  to 
o  o  o  o  o 

i^=40  ,  80  ,  and  140  ;  40  and  140  are  of  course  equidistant  from  the 
cylinder  axis. 

For  100  Hz,  Figure  IV. 10  shows  characteristic  dipole  patterns  for 

o 

all  three  elevations.  As  before,  >|)=80  provides  the  highest  levels.  The 
small  differences  between  the  \p=40°  and  the  140°  curves  suggest  that  for 
low  frequencies  the  rigid-body  structural  response  is  insensitive  to  the 
precise  axial  position  of  the  frame  being  driven. 

In  Figure  IV. 11  for  200  Hz,  the  departure  observed  from  the  dipole 
pattern  is  consistent  with  the  gap  between  the  two  curves  on  Figure  IV. 4 
widening  at  that  frequency.  The  1(1=80°  elevation  no  longer  dominates 
the  acoustic  field;  also,  the  differences  between  the  40°  and  140°  curves 
appear  more  pronounced  in  the  lower  hemisphere. 

.n  Figure  IV.  12  for  350  Hz  the  highest  levels  return  to  ijj=80°. 

However,  substantial  asymmetry  is  indicated  by  the  curves  for  the  other 
two  elevations. 

Figure  IV. 13  again  shows  a  dipole  pattern  at  all  elevations  with 
4>=80°  containing  the  highest  levels. 

Figure  IV.  14  shows  the  ip-  dependence  of  the  acoustic  field  for 
increasing  frequency;  the  500  Hz  curve  contains  the  largest  value  — 
roughly  -45  dB  at  i(i=100°;  350  Hz  follows  with  a  level  of  about  -51  dB 
at  ijj=60°  and  so  on.  The  predicted  on-axis  values  are  also  observed  to 

increase  with  frequency.  Interestingly,  these  values  seem  fairly  symmetric 

o  ... 

relative  to  t|j=90  even  if  the  rest  of  the  field  obviously  is  not;  i.e.  , 

two  far-field  observers,  one  facing  the  right  cap  and  the  other  the  left, 

would  hear  roughly  the  same  -69  dB  at  500  Hz, the  same  -78  dB  at  350  Hz,  etc. 
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F.  Iterated  Solutions:  Drive  on  Shell 
F. 1  Summary 

In  Chapter  III  iteration  algorithms  were  formulated  to 
solve  the  coupled  fluid-structure  systems  for  the  frame-stiffened,  finite 
cylindrical  shell.  It  will  be  recalled  that  for  n=0  the  resulting 
pressure  and  velocity  modal  vectors  are  equal  in  rank,  and  that  for  n^-1 
the  pressure  vector  is  two  entires  longer  than  the  velocity  vector  be¬ 
cause  cap  element  motion  is  then  taken  to  be  identically  zero.  For  both 
n=0,  n>l,  the  first  technique  was  the  standard  Jacobi  iteration  (J)  method 
adapted  to  coupled  systems.  As  previously  discussed,  Jacobi  iteration 
uses  a  single  iteration  counter  and  convergence  (or  divergence)  occurs 
after  applying  the  algorithm  a  prescribed  number  of  steps.  The  second 
form  of  iteration  for  both  n=0  and  n>l  equation  sets  is  actually  a  two- 
counter  process  which  uses  an  algorithm  nested  within  another.  As  was 
pointed  out  in  Chapter  III,  part  of  the  partitioned  subsystems  could  be 
solved  either  exactly  or  approximately  (Eqs.  III.C-8b  and  D-5c) .  Here 
we  use  the  exact  option  so  that  these  two-step  algorithms  fall  in  that 
category  of  linear-system  iterative  solution  techniques  which  Ref.  29 
terms  "semi-exact." 

Iterated  solutions  were  attempted  here  first  by  the  J  technique 
for  every  n  mode  of  every  frequency.  Successful  convergence  was  usually 
observed  for  the  lower  frequencies  between  20  and  100  Hz.  However,  as 
frequency  increased,  the  J  method  was  found  to  diverge  often.  For  such 
cases  the  two-step  algorithms  were  applied  with  both  counter  limits  set 
to  3  and  convergence  to  an  acceptable  solution  was  always  observed.  Thus, 
in  the  following  presentation  of  iterated  results  the  figure  label  "iteration 
counter (s) =3"  states  that  the  approximate  solution  was  obtained  either  in 
three  J  steps,  or  in  three  steps  of  the  two-step  processes;  "iteration 
counter(s)=l"  always  refers  to  the  solution  after  a  single  J  step,  which, 
as  discussed  in  Chapter  III,  is  the  quasi  in-vacuo  result  obtained  by 
approximating  the  structural  impedance  matrix  by  its  diagonal,  which  is 
then  inverted.  The  results  shown  were  chosen  as  representative  of  the 
body  of  data  generated. 
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F. 2  Iterated  Solution :  Frequency  Sweeps 

Figure  IV.  15  shows  exact  and  approximate  solutions  for  the 
frequency  sweep  at  i|»=80°,  4>=0  for  the  drive-on-shell.  The  exact  curve 
is  the  same  as  that  on  Figure  IV.  3  with  <j>=0.  Although  the  first  J  step 
yields  a  poor  estimate,  an  excellent  approximation  is  obtained  in  three 
generally  nested  steps  for  all  frequencies  except  perhaps  400  and  450  Hz. 

The  directivity  patterns  will  show,  however,  that  such  close  agreement 
was  not  obtained  at  other  orientations  in  the  acoustic  field. 

Figure  IV.  16  shows  the  frequency  sweep  at  the  two  far-field  positions 
coinciding  with  the  cylinder  axis  (<t>=0,ip=0  and  <j>=0,  ijj=180°)  .  Below  200  Hz 
even  the  one-step  solution  appears  adequate.  The  iterated  solution  using 
counter(s)  set  to  3  is  again  for  all  practical  purposes  the  exact  solution. 

F. 3  Directivity  Patterns:  Iterated  Solutions 

Figure  IV.  17  shows  exact  and  approximate  directivity  patterns 
on  the  plane  defined  by  ip=80°  with  <f>  varying  from  0  to  360°.  The  exact 
solution  is  taken  from  Figure  IV. 5  (different  scale).  We  note  that  the 
degree  of  improvement  vs  number  of  iterations  appears  to  be  global  in  that 
convergence  is  equally  good  for  all  values  of  circumferential  angle. 

Figure  IV.  18  shows  pattern  and  convergence  behavior  for  300  Hz.  The 
convergence  rate  for  this  case  appears  greater  than  that  for  100  Hz  since 
the  zeroth  solution  is  now  a  worse  first  estimate. 

In  Figure  IV.  19  for  400  Hz  convergence  ceases  to  be  globally  mono¬ 
tonic.  Thus,  after  three  steps  the  solution  for  <f>=180°  seems  to  be  down 
about  8  dB  and  that  for  $=0  up  4  dB;  this  4  dB  overshoot  was  seen  in 
Figure  IV.  15  at  400  Hz.  Figure  IV.  20  shows  the  pattern  for  500  Hz. 

Figures  IV.  21-24  show  directivity  patterns  on  the  plane  containing 

the  cylinder  and  <j)=0  axes.  For  all  figures  the  circumferential  angle  <j> 

o 

has  value  0  and  the  elevation  angle  ip  varies  from  0  to  180  .  Convergence 
was  observed  to  be  fairly  good  over  the  frequency  spectrum  investigated 
by  the  method. 
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It  should  be  pointed  out  that  within  the  frequency  range  investigated 
here,  i.e.,  below  500  Hz,  several  fluid-loaded  resonances  were  measured  in 
Ref.  40,  part  2;  and  that  although  only  far-field  predictions  have  been 
presented  here,  the  iterative  algorithms  were  observed  to  converge  for  all 
n  modes  to  the  correct  structural  response  (and  associated  near-field  sur¬ 
face  pressure)  within  this  resonant  range. 

G.  Comparison  to  Experiment 

Figure  IV.  25  shows  measured  beam-aspect  (4i=90°,  (J>=0)  levels 

for  the  cylindrical  structure  of  which  Figure  IV. la,  b  shows  one  of  five 
40 

equal  segments.  Also  shown  are  predictions  by  the  present  model  for 
frequencies  up  to  500  Hz,  and  by  a  simpler  infinite-cylinder  fluid  model 
for  frequencies  between  500  and  2500  Hz.  The  force-dipole  solution  is 
also  indicated. 

As  previously  pointed  out,  and  as  again  the  figure  shows,  predictions 
for  the  low  frequency  range  hover  about  the  force  dipole  with  more  con¬ 
sistent  departures  taking  place  as  frequency  increases.  Except  for  the 
measured  anti-resonance  dip  at  1  kHz  for  the  drive-on- frame  case,  pre¬ 
dictions  and  the  dipole  curve  would  seem  to  underestimate  the  mean  of  the 
data  by  about  25  dB. 

However,  the  results  of  previous  experimental  and  theoretical 
41  42 

studies  '  strongly  suggest  that  beam-aspect  acoustic  spectra  of  force- 

driven  cylindrical  structures  similar  to  that  investigated  here  are 

closely  approximated  by  the  force  dipole  at  the  lower  frequencies,  and, 

at  higher  frequencies  below  coincidence,  by  the  infinite  plate  solution 

42 

known  classically  to  be  6  dB  above  the  force  dipole.  Junger  for  example, 
has  identified  three  primary  frequency  ranges  for  a  typical  heavy-framed 
submarine  hull.  In  the  lower  of  these  the  ship  behaves  as  a  rigid  beam 
and  the  force-dipole  solution  is  roughly  valid.  For  the  next  higher, 
compartment-resonance,  range  he  has  established  structural  damping  criteria 
under  which  bounding  frames  play  a  subdominant  role  in  the  response  function 
which  then  corresponds  to  the  infinite-plate  solution;  although  at  fluid- 
loaded  resonances  peak  levels  may  otherwise  far  exceed  those  by  analytic 
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asymptotic  formulae,  he  finds  that  the  latter  nevertheless  provide 
good  mean-value  curves.  Finally,  above  the  hull's  ring  resonance  he 
determines  that  if  the  shell  plating  is  treated  as  an  infinite  system, 
its  acoustic  field  closely  approximates  the  frequency-averaged  measured 
data. 

Therefore,  the  discrepancy  in  overall  behavior  between  present 
and  earlier  measurements  plus  theoretical  conclusions  remains  to  be 
explained. 

H.  Computational  Times  for  Exact  vs  Iterated  Solutions 

We  end  the  discussion  here  with  a  brief  summary  of  the  time- 

savings  advantage  of  approximate  solutions.  In  all  fairness,  however, 

it  should  be  pointed  out  that  for  the  present  modal  formulation  total 

computational  costs  were  dominated  by  the  generation  of  matrix  elements 

for  the  fluid  system  rather  than  by  the  solution  technique  applied  to 

effect  its  inversion.  In  effect,  the  A  constant  in  Eq.  III. E-la  had  a 

large  value  as  a  result  of  the  integration  routines  used  here.  Neverthe- 

2 

less,  if  for  the  moment  the  AN  system-generation  contribution  is  ignored 
in  Eq.  III. E-la,'  one  may  calculate  part  of  the  N-x»  asymptotic  value  of  the 
savings  ratio  between  iterated  and  exact  solution,  i.e. ,the  (B/c)-l/N  part 
of  Eq.  III.B-lb  (N=M+2  here).  This  measure  of  savings  is  quantified  next. 

Table  IV. 3  shows  computational  times  needed  to  solve  a  given,  gen¬ 
erated  coupled  system.  Percent  time  savings  over  the  exact  solution  are 
tabulated.  The  first  of  the  approximate  solution  entries,  K^=l ,  K2=l , 
refers  to  the  J  method  with  and  as  defined  in  Chapter  III.  The 
second  entry,  K  =K  =3,  re^ers  to  t^le  nested,  two-step  iteration  process, 
for  which  a  set  of  K^,  combinations  were  also  applied  yielding  solutions 
of  accuracy  ranging  from  the  generally  acceptable  (K  '  to  *-he  9enerallY 

unacceptable  (Kj=K  =1)  . 

Had  a  technique  such  as  DAA  been  used  to  model  the  fluid  in  the 
present  problem,  the  ratio  of  total  matrix  inversion  to  matrix  generation 
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costs  would  have  been  essentially  given  by  (B/Ol/N.  The  reason  for  this 
of  course  is  that  DAA  requires  a  single  fluid  calculation  valid  for  all 
frequencies  and  thus  eliminates  much  of  the  system-generation  effort. 
However,  the  final  solution  would  then  be  doubly  approximate  because  the 
approximate  inversion  process  would  then  operate  on  an  approximate  fluid 
model. 

Finally,  as  discussed  in  Chapter  III,  SOR-type  recasting  of  fluid 
and  structural  systems  should  be  expected  to  improve  convergence  rates 
without  significant  computational  penalties.  Comparison  of  exact  and 
iterated  solutions  should  appear  then  even  more  in  favor  of  the  approximate 


method. 
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V.  SUMMARY  AND  CONCLUSIONS 

A  theoretical  model  has  been  developed  to  predict  acoustic 
radiation  from  a  submerged,  finite-cylindrical  shell  internally 
stiffened  by  ring  frames  and  with  flat,  movable  end  caps.  The  model 
has  been  applied  to  a  structure  for  which  recent  acoustic  data  has 
been  collected.  Predictions  differ  from  existing  measurements  by 
roughly  25  dB;  however,  they  appear  consistent  with  the  results  of 
a  number  of  previous  experimental  and  theoretical  studies. 

Iterative  schemes  have  been  developed  to  solve  the  associated 
fluid-structure  coupled  systems  over  a  moderate  range  of  frequencies 
encompassing  several  fluid-loaded  resonances.  We  have  demonstrated 
that  such  techniques  may  be  applied  to  generate  essentially  the 
exact  solution  in  a  reasonably  small  number  of  iterations.  By  com¬ 
paring  computational  costs  to  those  for  the  exact  result,  we  found 
iteration  to  be  economically  advantageous  for  all  frequencies  considered. 

Lastly,  we  have  noted  that,  for  an  arbitrary  fluid-loaded  struc¬ 
ture  modeled  using  finite  differences  or  elements,  the  ratio  of  total 
computing  cost  by  iteration  to  that  using  exact  inversion  should  de¬ 
crease  with  system  size.  The  development  and  implementation  of  tech¬ 
niques  such  as  those  used  here  is  therefore  recommended. 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC, 


(1) 

Tt 

/ 

0 

where 


where 


APPENDIX  II. A 


2 

a  I 


2  (-1) 


2n 


ri  2  ,  -.2.3/2 

{1+r  +(z-z)  } 


dacosna  (1+rcosa)  [1-ika  /l+x^+  (z-z)  2  /l+£^cosa]  exp{ika  Jl+ r  +(z-z)  /l+^^cosa} 


{l+I^cosa}3^2 


;  =  2r/[l+r2+(z-z)2] 


al 


2  (-1) 


n  it 


3n 


/ 


/~2  CS  2*  / - - - 

dacosna  exp{  ika  /r  +r“+ ( z+L/2a)  /!+£  cosa} 


i^:2+r2+  ( z+L/2a)  2  0 


/l+C2cosa 


a3I 


2  (-1) 


4n  {r2+r2+(z+L/2a)2}3/2 


J*  dacosna 
0 


[l-ika/42+r2+  (z+L/2a')  2  /l+E  cosa]  exp  { ikav42+r2+  (z+L/2a)2 A+£. 


oosa} 


{ l+^cosa} 


3/2 


-2-2  2 
=  2rr/[r  +r  +(z+L/2a)  ] 


2  (-1) 


//  2  n  2 

r  +r  +(z-L/a)  /I+]f~~cosa} 


al 


3n  /r^+r2+ (z-L/2a) 2 


/ 


/T+^jCosa 


0 


PLANE  NORMAL  TO  SHELL  AXIS 
DRIVE  ON  SHELL.  F=200  Hz 


r 4  a 

///' 

/  <T  /  ap 


Ot  \  \o 

TO 

\  X  A  > 


nA.  \ 


-Ml'' 


LD  LD  LD 
Ll)  U)  Ul 

o  n  n 

□  □  □ 

t\)  ^  00 


LD  LD  LD 
Z  Z  Z 
<C  <  < 


>  >  > 

Lxl  Ul  Ul 
_l  _J  _) 
Ul  UJ  Li] 


Figure  IV. 6  Predicted  Acoustic  Directivity  vs  <p  for  200  Hz  with  Point  Drive  on  Shell  Middle 


-70 


Figure  IV.  3  Calculated  Spectra  of  Radiated  Sound  at  (<j>=0,  ^=80  )  ,(<{>=180  ,^=80  )  for  Drive  on  Shell 
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TABLE  IV. 2 

STRUCTURAL  AND  ACOUSTIC  WAVENUMBERS,  AND  REQUIRED  SYSTEM  SIZE 


f  (Hz) 

kfa(s,188/f) 

M *k  a) 

7T  f 

ka(=f/1000) 

20 

.84 

6 

.02 

50 

1.  32 

9.3 

.05 

100 

1.88 

13.2 

.1 

200 

2.66 

18.6 

.2 

500 

4.20 

29.4 

.5 

1000 

5.95 

41.6 

1. 

2000 

8.41 

59.0 

2. 

5000 

13.  30 

93.0 

5. 
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TABLE  IV. 1 

CYLINDER  STRUCTURAL  PARAMETERS 


Cylinder  length 
Shell  thickness 
Shell  radius 
Caps'  thickness 
No.  of  frames 
Inter-frame  spacing 
Frame  height 
Frame  width 


230" 

.3125"  (Aluminum) 
10. 344" 

1.5"  (Steel) 

4  (Aluminum) 

46" 

1.0625" 

3.25" 


Assumed  Loss  factor  for  all  structures 


.05 


A*  /Kitw'nu**  j/^/f 


Idealized  Geometry  Used ' for  Theoretical  Predictions,  Structural  Parameters  are  Indicated 
for  Shell,  Caps,  and  Frames  (also  Table  IV. 1)  Drives  considered  here  are  shown:  on  shell 
on  frame. 


TOP  VIEW 


ACCESS  1—  7"  8"  — |—  8"  — |—  7" 

MOLE 


Figure  IV.la,b  From  Ref.  40;  Top  and  Side  Views  of  a  Typical  AL  Interframe 

Cylinder  Section  of  the  Experimental  Model  (AL=46"  in  Fig.  IV. 2) 


67 


Figure  II. 1  Finite  Cylinder  Model  Showing  (a)  Side  View  with  (r,z)  Origin  and  Structural 

Dimensions  for  Shell, Caps,  and  Frames;  (b)  Coordinate  System  for  Far-field  Calculations 
(c)  Cylinder  Cross-section  Showing  Structural  Ring-frame  Parameters. 
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where  "det"  denotes  the  determinant  of  the  matrix  in  Eq.  II.B-4. 
(2)  Modal  Impedances  for  Ring  Frames 


~  28 
Modal  impedances  Z  f  etc.,  have  been  given  Harari  in  dimensional 
nil 


form  and  relate  net  frame  forces  to  frame  displacements  rather  than  velocities. 


,H 


The  relationship  between  Harari' s  "impedance"  Z  and  our  normalized  quantity 

nil 
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nil 
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nil  ka  2 
pc 


and  similarly  for  the  rest.  Thus  they  are  given  by 

4  2  2. 
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n31 


fi. 


bS 


in  [l  -  (l  - 1^-)  4  <w2] 


where  b  is  the  frame  width  as  defined  in  Fig.  II.  2,  which  also  shows 
is  the  frame  material  density,  c_  its  compressional  wave  speed,  and  its 


structural  loss  factor;  also, 


fi,  =  (wa/c  )//l-in. 


In  the  above  equations  all  spatial  quantities  appear  normalized  by  the  shell 
radius  a. 
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,1+v,  /iniTa, 
A2  =  (— )  n  (— > 


a  .nvrra. 

A3  =  v  (— ) 


A4  =  A2 

a5  -  n2-a+e2){n2  +  (i^)(=i)2} 


r .  2  .  2  . mir 9,2. 
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where 


n  =  (ua/c  )  /  A- in  ; 
P 


c  =  shell  compressional  wave  speed; 
P 


=  structural  loss  factor;  and 


6  =  (h/a)/AT 


The  constants  By  By  By  Bg  are  given  by 


WVs 


A2A6~A3A5 
3  det 


R  vw 

B7  det 


A. A  -A. A . 
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APPENDIX  II. B 


(1)  Structural  Operators  and  Constants 

The  thin-shell  theory  operators  L  ,  ...,L  in  Eqs.  II.B-la  are  given 

1  9 

by  (Refs.  34,  and  12). 


L  =  _ll  +  A ZV  il.  + 
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The  constants  A  , ...»  A  in  Eq.  II.B-4  are, 
1  9 


.  2  .imra,  ,1-v.  2 

A1  =  n  -  (— )  -  (— )n 
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Figure  IV.  9  Predicted  Acoustic  Directivity  vs  i|>  for  <j>=0;  Point  Drive  on  Shell  Middle. 
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Figure  IV.  13  Predicted  Acoustic  Directivity  vs  ip  for  500  Hz  with  Point  Drive  on  Near-middle  Frame;  ip=40  ,  80  ,  140 
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Predicted  Acoustic  Directivity  vs  tp  for  <f>=0;  Point  Drive  on  Near-middle  Frame 
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Figure  IV.  23  Comparison  of  Exact,  and  Approximate  Directivity  Patterns  using  Iteration;  <t>=0,  ij>  variable; 
f=400  Hz;  Drive  on  Shell 
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TABLE  IV.  3 

COMPUTATIONAL  TIME  COMPARISON  BETWEEN  EXACT  AND  ITERATED  SOLUTIONS  VS. 


ITERATION  COUNTERS 


V  *2 


Solution  #  K2 


Execution  Time 
(secs) 


%  Time  Savings 
Over  Exact 


SOUND  PRESSURE  LEVEL-tlB  RE  l/i  Pa  AT 
lm  PER  NEWTON  INPUT  FORCE 


> 


I 


\ 
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Figure  IV.  25  Comparison  of  Predicted  to  Measured  Spectrum  (from  Ref.  40)  at 
(j>=0,  <J)=90°;  Predictions  for  f>500  Hz  use  an  Infinite  Cylinder 
Model  for  the  Fluid.  94 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC. 


REFERENCES 


1  R.  Martinez  and  J.  M.  Garrelick,  "An  Exact  Analysis  of  the  Force- 
Excited  Submerged  Spherical  Shell  using  Iterative  Techniques," 

CAA  Rept.  C-911-288,  Feb.  1983. 

2  L.G.  Copley,  "Fundamental  Results  Concerning  Integral  Represen¬ 
tations  in  Acoustic  Radiation,"  JASA  44  (1) ,  1968,  pp.  28-32. 

3  H.A.  Schenck,  "Improved  Integral  Formulation  for  Acoustic  Radiation 
Problems,"  JASA  44  (1),  1968,  pp.  41-58. 

4  K.  Brod,  "On  the  Uniqueness  of  Solution  for  all  Wavenumbers  in 
Acoustic  Radiation,"  76  (4),  1984,  pp.  1238-1243. 

5  I.  Mathews  and  D.  Hitchings,  "Acoustic  Radiation  from  Three- 
Dimensional  Elastic  Structures,"  Innovative  Numerical  Analysis 
for  the  Engineering  Sciences,  UVA  Press,  1980,  pp.  69-78. 

6  G.T.  Schuster  and  L.C.  Smith,  "A  Comparison  Among  Four  Direct 
Boundary  Integral  Methods,"  JASA,  77  (3),  1985,  pp.  850-864. 

7  H.  Huang,  "Helmholtz  Integral  Equations  for  Fluid-Structure 
Interaction,"  Advances  in  Fluid-Structure  Interaction  -  1984, 
presented  at  the  1984  Pressure  Vessel  and  Piping  Conference  and 
Exhibition,  San  Antonio,  Tx. ,  June  17-21,  1984. 

8  G.  Chertock,  "Sound  Radiation  from  Vibrating  Surfaces,"  JASA 
36  (7),  1964,  pp.  1305-1313. 

9  R.  Courant  and  D.  Hilbert,  Methods  of  Mathematical  Physics, 

Vol.  I,  Interscience,  1966,  p.  136. 

10  B.E.  Sandman,  "Numerical  Fluid  Loading  Coefficients  for  the  Modal 

Velocities  of  a  Cylindrical  Shell,"  Comp.  &  Struc. ,  6,  1976,  pp. 
467-473. 


9 


i 


95 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES.  INC. 


11  B.E.  Sandman,  "Fluid-Loading  Influence  Coefficients  for  a  Finite 
Cylindrical  Shell,"  JASA  60  (6),  1976,  pp.  1256-1264. 

12  M.C.  Junger  and  D.  Feit,  Sound,  Structures,  and  Their  Interaction, 
1st,  ed. ,  MIT  Press,  1972,  p.  207. 

13  V. A.  Kozyrev  and  E.L.  Shenderov,  "Radiation  Impedance  of  a  Cylinder 
of  Finite  Height,  "Sov.  Phys.  Acous.,  26  (3),  1980,  pp.  230-236. 

14  V.  A.  Kozyrev  and  E.L.  Shenderov,  "Radiation  Impedance  of  a  Cylinder 
with  Acoustically  Compliant  Ends,"  Sov.  Phys.  Acous.,  27  (3),  1981, 
pp.  247-248. 

15  J.S.  Patel,  "Radiation  and  Scattering  from  an  Arbitrary  Elastic 
Structure  using  Consistent  Fluid  Structure  Formulation,”  Computers 
&  Struct.,  9,  1978,  pp.  287-291. 

16  T.L.  Geers,  "Transient  Response  Analysis  of  Submerged  Structures," 
Finite-Element  Analysis  of  Transient  Nonlinear  Structural  Behavior, 
edited  by  T.  Belytschko,  et.a].,  ASME,  N.Y.,  1975,  pp.  59-84. 

17  T.L.  Geers,  "Doubly  Asymptotic  Approximations  for  Transient  Motion 
of  Submerged  Structures,"  JASA,  64(5),  1978,  pp.  1500-1508. 

18  T.L.  Geers  and  C.A.  Felippa,  "Doubly  Asymptotic  Approximations  for 
Vibration  Analysis  of  Submerged  Structures,"  Lockheed  Palo  Alto 
Research  Laboratory  Rept.  prepared  for  DTNSRDC/Annapolis ,  Sept.  1981. 

19  H.  Huang  and  Y.F.  Wang,  "Asymptotic  Fluid-Structure  Interaction 
Theories  for  Acoustic  Radiation  Prediction,"  Advances  in  Fluid-*- 
Structure  Interactions  -■  1984,  see.  ref.  7  for  details. 

20  D.  Ranlet,  F.L.  Dimmagio,  H.H.  Bleich  and  M.L.  Baron,  "Elastic 
Response  of  Submerged  Shells  with  Internally  Attached  Structures 
to  Shock  Loading,"  Comp.  &  Struct.,  7,  1977,  pp.  355-364. 


96 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC, 


21  L.H.  Chen  and  D.  G.  Schweikert,  "Sound  Radiation  from  an  Arbitrary 
Body,"  JASA,  35  (1),  1963,  pp.  1626-1632. 

22  M.L.  Baron,  L.A.  Matthews,  and  H.  H.  Bleich,  "Forced  Vibrations  of 
an  Elastic  Circular  Cylindrical  Body  of  Finite  Length  Submerged  in 
an  Acoustic  Medium,"  Weidlinger  Tech.  Rept.  1,  prepared  for  ONR, 
June  1962. 

23  J.  M.  McCormick  and  M.L.  Baron,  "Sound  Radiation  from  Submerged 
Cylindrical  Shells  of  Finite  Length,"  Weidlinger  Tech.  Rept.  3, 
prepared  for  ONR,  June  1964. 

24  J.  M.  McCormick  and  M.  L.  Baron,  "Sound  Radiation  from  Submerged 
Stiffened  Cylindrical  Shells  of  Finite  Length  with  Flexible  End 
Caps,"  Weidlinger  Tech.  Rept.  9,  prepared  for  ONR,  Aug.  1970. 

25  D.  G.  Crighton,  "Resonant  Oscillations  of  Fluid-Loaded  Struts,” 

J.  of  Sound  &  Vib.,  87  (3),  1983,  pp.  429-437. 

26  C.B.  Burroughs,  "Acoustic  Radiation  from  Fluid-Loaded  Infinite 
Circular  Cylinders  with  Periodic  Ring  Supports,"  JASA,  75,  (3), 
1984,  pp.  715-722. 

27  J.E.  Cole,  "Parametric  Study  and  Design  Implications  of  the  Low 
Frequency  Response  and  Radiated  Noise  of  a  Submarine  Hull,"  CAA 
Rept.  C-938-274. 13E ,  prepared  for  Naval  Sea  Systems  Command, 

Aug.  1983. 

28  A.  Harari ,  "Wave  Propagation  in  Cylindrical  Shells  with  Finite 
Regions  of  Structural  Discontinuity,"  JASA,  62  (5) ,  1977,  pp.  1196- 
1205. 

29  D.M.  Young,  Iterative  Solution  of  Large  Linear  Systems,  Academic 
Press,  N. Y. ,  1971. 

30  M.K.  Au-Yang,  "Natural  Frequencies  of  Cylindrical  Shells  and  Panels 
in  Vacuum  and  in  a  Fluid,"  JASA,  57  (3),  1978,  pp.  341-355. 


97 


CAMBRIDGE  ACOUSTICAL  ASSOCIATES,  INC 


Y.W.  Chang  and  C.Y.  Wang,  "An  Eulerian  Method  for  Large  Displacement 
Fluid  Structure  Interaction  in  Reactor  Containment,"  Computational 
Methods  for  Fluid-Structure  Interaction  Problems,  AMD  Vol.  26, 

Winter  Annual  Meeting  of  the  ASME,  Atlanta,  Ga.,  Nov.  27-Dec.  2, 

1977,  pp.  1-14. 

T.  Belytschko,  "Fluid-Structure  Interaction,"  Comp.  &  Struc. ,  12, 

1980,  pp.  459-469. 

W.  Whitlow  and  W.  L.  Harris,  "Lifting  Airfoils  in  Transonic  Flow 
in  Parametric  Differentiation,"  AIAA  Paper  80-1394.  Presented  at 
the  13th  AIAA  Fluid  and  Plasma  Conf . ,  Snowmass  Col.,  July  14-16,  1980. 

O.D.  Kellog,  Foundations  of  Potential  Theory,  Dover,  N.Y.,  1953,  p.  299. 

A.  W.  Leissa,  Vibration  of  Shells,  NASA  SP-288,  1973#  for  Donnell- 
Mushtari  theory  see  p.  32. 

M.E.  Goldstein,  Aeroacoustics ,  McGraw-Hill,  N.Y.,  1976,  p.  119. 

A. A.  Cannarozzi  and  E.  D'Anna,  "A  Routine  for  Solving  Sparse  Linear 
Systems  in  Structural  Analysis,"  Comp.  &  Struc.,  6,  1976,  pp.  99-105. 

H.  Peterson  and  E.  P.  Popov,  "Substructuring  and  Equation  System 
Solutions  in  Finite-Element  Analysis,"  Comp.  &  Struc.,  7,  1977, 
pp.  197-206. 

A.  Maher,  P.G.  Kessel,  and  R.D.  Cook,  "A  Partitioned  Finite-Element 
Method  for  Dynamical  Systems,"  Comp.  &  Struc.,  18,  pp.  81-91. 

C.  B.  Burroughs,  S.I.  Hayek,  D.A.  Bostain,  and  J.E.  Hallander, 

"Acoustic  Radiation  from  Single  and  Double  Ribbed  Circular  Cylin¬ 
drical  Shells,"  Penn.  State  ARL  TM  84-76,  March  1984;  also,  "Resonant 
Frequencies  and  Mode  Shapes  for  Single  and  Double  Cylindrical  Shells," 
Penn.  State  ARL  TM  83-169,  Sept.  1983. 

H.B.  Ali  and  D.  Feit,  "Acoustic  Measurements  and  Prediction  for  a 
Submarine  Model  -  A  Heuristic  Appraisal,"  DTNSRDC  Rept.  78/042, 

May  1978. 

M.C.  Junger,  "Force  Dipole  Noise  Radiation  by  Submarine  Hulls," 

CAA  Rept.  C-890-274,  13A  prepared  for  Naval  Sea  Systems  Command, 

July  1983. 


END 

FILMED 

7-85 


DTIC 


